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ABSTRACT 


The objective of this work has been to investigate the feasibility 
of microwave holographic Imaging of targets near the earth usxng a 
large random conformal array on the earth's surface, the array and target 
volume being illuminated by a CW source on a geostationary satellite. 

Specific goals are defined in section i. Introduction, where a summary 
of the report is presented and the main conclusions are outlined. High- 
lights of the work accomplished are as follows. Section 2 establishes the 
necessary geometrical formulation for illuminator- target-array relation- 
ship calculations. This is applied in section 3 to the calculation of sig- 
nal level3 resulting from L-band illumination supplied by a satellite 
similar in capability to the ATS-6. The need for extensive spatial and 
temporal processing capability (as would be provided with multiple an- 
tenna elements together with multiple narrow-band filtration of signals) 
to build up sufficient signal-to-noise ratio of signals reflected from 
space and airborne vehicles is indicated. It is found that for the para- 
meter ranges assumed, a system for surveillance of airborne targets could 
result in a practically feasible system. 

In section 4 the relations between direct and reflected signals are 
analyzed and the composite resultant signal seen at each antenna element 
is described. Processing techniques for developing directional beam for- 
mation as well as SNR enhancement are developed in section 5. Spatial 
processing (combination of signals from the individual antenna elements) 
techniques based on either coherent detection or square-law detection, are 
shown possible and comparative evaluations are presented. Temporal pro- 
cessing techniques appropriate to an aircraft target environment are ob- 
tained. Also presented in this section are a representative system block 
diagram and a discussion of possible nulling antenna applications. 

Section 6 develops the angular resolution and focusing characteris- 
tics of a large array covering an approximately circular area on the 
ground. In section 7, the necessary relations are developed between 
the SNR achievable by coherent integration subsequent to beam formation, 
the range, speed, and radial acceleration of the aircraft targets, and 
the size and number of elements in the array. 

Numerical results are presented in section 8 for a possible air 
traffic surveillance system. It is shown that a system with several 
thousand elements spread over an area several hundred wavelengths in dia- 
meter should be feasible. Representative antenna element characteristics 
are defined-, and the utility of a pulse-radar to aid in focusing as well 
as provide explicit target range information is described. 

Finally, In section 9, a simple phase correlation experiment is de- 
fined that can establish how large an array may be constructed in the 
L-band of interest working in conjunction with the ATS-6 satellite. It 
is anticipated that carrying oat this experiment will also be useful in 
developing techniques for design and evaluation of large spaceborne antennas. 
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RESEARCH IN LARGE ADAPTIVE 
ANTENNA ARRAYS 

1. INTRODUCTION 

1.1 Objectives and Goals 

The objective of this work has been to Investigate the feasibility 
of microwave holographic imaging of targets near the earth using a large, 
random, conformal array on the earth's surface, the array and the target 
volume being illuminated by a CW source on a geostationary satellite. 
Specific goals to be pursued are as follows: 

a. Development of system concept involving use of satellite signal 
as phase reference for array and also to illuminate target space. 

b. Determination of signal level requirements taking into account 
satellite transmitter power and gain, and forward scatter cross- 
section characteristics of typical aircraft and low-orbiting 
space vehicles. 

c. Formulation of algorithms for accurate angular location of tar- 
gets by coherent processing of received signals, accounting for 
Doppler displacement of scattered signals as well as angular 
and range location of targets. 

d. Determination of applicability of (possibly adaptive) null form- 
ing techniques in the array for reducing the apparent level of 
the direct signal from the satellite so as to minimize the dynamic 
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range between the satellite signal and signals scattered by tar- 
gets . 

e. Investigate the system utility that would be implied by the use 
of a separate ground-based pulse-radar to provide range and 
coarse angle information for targets of Interest. 

f. Determine significant frequency band tradeoffs comparing L-, 

S-, C-, and X-band, or higher, operating frequencies . 

g. Determine the feasibility of phase synchronization of a large, 
random, conformal array by signals from a geostationary satellite, 
developing a plan for an experiment that can lead to such de- 
termination. 

1.2 Summary 

In section 2 to follow the basic satellite - array geometry is de- 
scribed and the necessary formulas for calculating observation distances 
and angles are presented. 

Next (section 3) are covered the signal levels which result for both 
direct satellite - array and indirect satellite-target-array propagation 
paths ; assuming ATS-6 in geosynchronous orbit as the L-band illuminating 
source, typical airborne or low-orbiting space vehicles as targets, and 
antenna array elements located in the vicinity of our Valley Forge Re- 
search Center. A calculation is made of the number of array elements re- 
quired to obtain a detectable reflected signal; it is concluded that if 
no more than a few thousand elements can be used, it is necessary to limit 
the application to the surveillance of airborne targets ( not feasible for 
detecting space vehicles unless orders of magnitude more illuminating 
power were available) . 
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In section 4 explicit expressions for direct and reflected signals 
are presented so as to obtain a mathematical representation of resultant 
waveforms at each antenna element. 

Section 5 is concerned with the processing of signals, spatial and 
temporal combination as required for beam forming, which is the equiva- 
lent of microwave holographic Imaging of the illuminated volumes. Both 
coherent detection and incoherent square-law detection techniques are de- 
fined; it is concluded that square- law detection would be much simpler to 
Implement although it would suffer 3 decibels SNR degradation relative to 
the coherent detection schemes. A system block diagram is proposed, sug- 
gesting the feasibility of performing temporal filtering (FFT or doppler 
filter banks) subsequent to multiple beam formation. Finally a technique 
is presented for nulling out the direct ray so as to permit more effective 
detection of target echoes. 

Section 6 is concerned with the formation of beams by the large cir- 
cular array, including angular resolution characteristics as well as 
range-dependent focusing requirements. It is concluded that a coarse 
estimate of target range will be sufficient to permit near optimum angular 
resolution. 

in section 7 the dependence of signal detection capability on target 
range and speed and on spatial and temporal signal process i.ig parameters 
is obtained. Graphs arc obtained, presenting parameter dependences, in- 
cluding the effects of rate of change of range rate during observation in- 
tervals and showing dependence on carrier wavelength. 

In section 8 specific sets of numerical values are given, illustrat- 
ing possible system performance capability. 


tssffW.W 
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Pinally, in section 9 an experiment is described that can be considered 
as a first step in the implementation of the concepts developed. 

1.3 Conclusions 

The goals enumerated in section 1.1 above have been essentially car- 
ried out as planned. In general, the concept of microwave holographic im- 
aging of airborne vehicles for air traffic surveillance has been shown to 
be technically feasible. If it were decided to proceed with the develop- 
ment of the system, the next step would be to perform a more detailed sys- 
tem design as well as more detailed comparison with other techniques for 
obtaining the deal, id target information. 

It is expected that adaptive nulling techniques can be usefully 
applied; also, the experiment planned for the next year can provide addi- 
tional environmental information as to the applicability of the concepts 
developed. 

2. SATELLITE- ARRAY GEOMETRY 

As a convenient point of departure, we have decided to employ selected 
characteristics of the ATS-6 satellite and its systems [1,2] as a starting 
point in our system development. This will serve to block out regions of 
feasibility of the concepts developed and can lead to a reasonable indica- 
tion of hardware requirements. 

Consider a satellite in synchronous orbit radiating a 1550 MHz signal 


[1] The ATS-F&G Data Book, Goddard Space Flight Center, Greenbelt, Md., 
Revised edition, Sepatember 1972. 

[2] X-460-74-232 Applications Technology Satellite ATS-6 In Orbit Check- 
out Report, August, 1974, Goddard Space Flight Center, Greenbelt, Md. 
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(CW sinusoid unmodulated) from an antenna 30 feet ia diameter (parabo • 

loidal reflector) pointing at a random array of antenna elements on the 

earth's surface (for example at our Valley Forge Research Center site). 

Since the wavelength is some 20 cm or about 0.66 feet, the beamwidth of 

the 30 foot antenna will be about .022 radian*. Thus in the vicinity of 

the earch's surface, a volume in space will be illuminated that will be 

roughly cylindrical, concentric with the line of sight from satellite to 

antenna, having a diameter of .022 x 20,000 or about 440 nautical miles. 

The parameters of the reflected signals measured at a receiving array 

element (amplitude, frequency and phase) are functions of the distances 

involved (source-target distance and target-array distance) and their 

rates of change. Therefore it will be useful to find some expressions 

concerning the satellite-target-array geometry. 

. Let us denote, as shown in Figure 1, the distance between t 

. stationary satellite and the center of the earth by R , the radius of the 

6 

earth by and the distance between the satellite and a reference point 

R within the array by d • Let $ and ip be the longitude and the latitude 

£8 £ £ 

of the reference point R,and <p the longitude of the satellite position. 

s 

(The latitude of a geostationary satellite is ip m 0.) We are interested 

3 

in the distance between the satellite and the point R (i.e., d ) and the 

sr 

direction of the satellite with respect to the point R, i.e., the com- 
ponents of a unit vector k pointing towards the satellite in a suitably 
chosen rectangular coordinate system (see the system XYZ-R in Figure 1.) 

The distance d g ^ is easily found by 

d ■ £ (x - x )^ + (y - y )^ + (z - z )^]** 
sr 1 os or' w os 'or' v os or' J 


( 1 ) 
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where 


FIGURE 1. SATELLITE-ARRAY GEOMETRY 

*or a R r cos\(/ r C08(|» r 

y Qr - R r costJ» r 8in<^ r 


z ■ R sirvto 
or r r r 


X ■ R cosil> co&6 *= R cos<fc 

08 8 a y 8 8 8 


y 0 « “ V 08 ^ sln K " V in *s 


Z 08 " R 8 Sin ^S " 0 


( 2 ) 


are the coordinates of R and S in the rectangular coordinate system 

x y z -0, with z -axis passing through the North pole of the earth and 
0 0 0 0 


the x Q -axis passing through Greenwich meridian (see Figure 1) 


IlLVliv 'K I, '' ! .IT’i t ! i' I’ 
ORIGCNAL i'.' : • !><*.>*. 
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The components k , k and k of the vector k in the XYZ-R coordi- 
xr yr zr 

nate system, where the X-axis points towards south, the Y-axis points 
towards east and the Z-sxis is normal to the earth surface, are equal re- 
spectively to the components k fl , k, and k in the spherical coordinate 

o $ r 

system 6$R-0. Therefore 


k 


xr 


k cos9 cos* -I- k cos6 sin* - k sinO 
xo r T r yo r T r zo r 


k _ ■ k . “ -k sin$ + k cos$ 
yr $ xo r yo r 


k “ k ■ k sin9 cos* + k sin9 sin* + k cos0 

zr r xo r T r yo r r zo 


(3) 


Here, k , k and k are components of the vector k in the X Y Z -0 co- 

xo’ yo zo ^ o o o 

ordinate system, i.e.. 


k 


xo 


X - x 
os or 

d sr 


k 


yo 


y - y 

^os ^or 


sr 


k 


2 -Z "2 

os or or 

cr - — — - „ rx ■ 

zo d d 

sr sr 


(4) 


and 


9 f ■ ^ (as can be seen in Figure 1). 

* 

For example if the satellite is somewhere above the Galapagos Islands 
(say * » -94,5°) and the array is in the Valley Forge area ($ * -75.5° 

8 T 

and <J» - 40°), by using (1), (2), (3), and (4) we get d . - 37.8106 x 10 3 km, 

k a 0.6779, k = -Q.3629 and k = 
xr yr zr 


0.6394. 



Hence, the elevation angle a and the azimuthal angle 4 of the satellite 

w ill be given by: a ■ v/2 - arcsing ■ 39.747 degrees and 

si r 

4 ■ arctan(k /k ) * -28.161 degrees (i.e., 28.161 degrees south-west), 
yr xr 

3. SIGNAL LEVEL CALCULATIONS 

Consider the signals received by a typical antenna element. First 
there will be a sinusoid received as a plane wave from the direction of 
the satellite. Next there will be a sinusoidal echo from each target 
present in the illuminated volume. The frequency of the echo will be dis- 
placed from that emitted by the satellite by the Doppler shift given by 
v 

— f where v is the rate of change of propagation pathlength (satellite 
w o r 

to target to antenna). For f ™ 1550 MHz, this becomes f. B 10.3 v Hertz 

o a r 

where v f is in m/s. Typical v^ values are 1000 m/s for fast moving air- 
craft, 7500 m/s for low orbital satellites. Corresponding f^ values are 
thus 10.3 KHz and 77 KHz. Such doppler shifts are important parameters 
that must be considered in any schemes for signal enhancement or informa- 
tion extraction. 

3.1 Preliminary Environmental Considerations 
A. Formulas required for signal level calculations 
For the direct link [3] we have: 

1. (S/N)* (P ) + (G ) + (Grd) + 2(5.)- 2 (R) - (B) - (NF q ) 
and for the reflected link (4]. 

2. (S/N) r - (P fc ) + (G t > + (Grr) + 2(X) + (a) 

- 2(R t ) - 2(R r ) - (B) (NF o > - (L) 


- (L) + 77 


(5) 

( 6 ) 


[3] R. Berkowits (cd.), Modern Radar . John Wiley & Sons, Inc., p. 13, 1965. 

[4] Ibid., p. 12. 
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where 

(P fc ) » transmitted power in dbw 

(G t ) ■» gain of transmitting antenna in db 

(Grr),(G ) “ gain of receiving antenna in db 

rd 

(X) ■ wavelength in db cm 

(&) ® distance of satellite to receiver in db nmi 

(X £ ) * distance of satellite to target in db nmi 

(R ff ) ■ distance of target to receiver in db csi 

(a) “ bistatic radar cross section of target in dbm^ 

(B) a signal bandwidth in db Hz 

a t ) ° loss factor in db 

(NF) ■ noise factor in db 

o 


B. Initial selection of representative numerical values 

Parameters similar to those of the ATS-6 satellite will be assumed. 
Basically, we have 40 watts transmitted at L-band (1,550 MHz) over a 
30 foot parabolic reflector antenna. The transmitter is assumed in syn- 
chronous orbit pointing at the receiver array assumed to be in the Phila- 
delphia- Valley Forge area. We can calculate as follows for each parameter: 


<P t > 

<G t ) 

(X) 


(40) - 116.02 dbw 


38.5 db| (ATS-6 book) 

,3x30 x100, , 

( q) - (19.35 cm) 

1.55x10* 


(L) ■ 4 db (typical value) 


(NF ) • 5 db 1 (typical value) 

o l - ■ 

(R) ■ 43.10 dbntni (see calculations Figure 1) 

(R ) *= (R) typical assumption assuming targets in illuminated 

region 





10 - 


C8r) 

(Grr) 

(Grd) 

(<0 


(100 nai) 
[ SdB^j 


]20 dbnail 


10 db® (nominal asswpeloo typical for airborne 
and apace vehicle sises) 

CL.6 x 10 5 Bs) » (52.04 dbHzt 


C. Single element S/B calculations 

For the direct link we calculate from equation (5) as follows: 


<p t ) 

- 16.02 

2(1) - 

86.20 

«£> 

- 33.50 

_JB> - 

52.04 

(CrI) 

» 5.00 

(UP ) - 

5.00 

2(A) 

» 25.74 

(2) ° 

4. Op 


162.26 
- 147.24 


147.24 


|l5.02dB 

- (S/N) d 



For the reflected link we calculate fro® (6) as follows: 


<V 

a 

16.02 

2(R t > 

- 

86.20 


m 

38.50 

20^) 

- 

40.00 

(Grr) 

■ 

5.00 

(B) 

a 

52.04. 

2(A) 

a 

25.74 

(HP„) 

a 

5.00 

(o) 


10.00 

a) 

4 a 

4.00 



95.26 



187.24 


- 

187.24 





E 

91.98 dB - 

(S/S) 




Note the apparent 107dB difference between direct and reflected rays. 
This is given more accurately by 


3- S|»| d 

Spffr 


^ ,R r 5 rd 


106.34 dB 


(Note here Is in m» not nmi) 


Rff RODUCffiiLiTYorraE 

ORIGINAL PaGE IS POOR 




The conclusion from this discussion is that while single** element S/N 
ratios will be satisfactory for direct link signal detection, the reflected 
l in k signal necessitates a combination of higher transmitted power, spatial 
processing (coherent combination, or its equivalent of signals from multiple 
antenna elements "beam forming") and temporal processing (time averaging) 
accomplished after doppler filtering. 

3.2 MINIMAL SIZE AND EFFECTIVE AREA OF THE ARRAY APERTURE 

Imaging of the volume of interest will be performed by effectively 
forming, from the same sets of array element data, many beams simultan- 
eously. Beam spacing will be close enough so that a target is within the 
3 dB width of one of them or, as is usually said, within a "resolution 
cell" . 

Presence of a target in a resolution cell will be assessed on the 
basis of the received energy reflected from the target during its stay 
in that resolution cell. The larger the received energy, the higher the 
probability that error will not occur. 

The amount of received energy depends on the transmitted power, the 
target characteristics, the effective receiving area and the required 
angular resolution. 

Limited illuminating power restricts the domain of applicability and 
the performance of our system. In this section we find simple relations from 
which, given the transmitted power, one can determine the required number 
of array elements so that certain targets can be handled and, conversely, 
the limitations on the target characteristics for a given number of array 
elements and given angular resolution. We also examine the effect of the 
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radlatloa wavelength on the required number of elements under the assump- 
tion of constant transmitting antenna gain. 

A fundamental quantity In our discussion will be the "effective re- 
ceiving area" or "effective area of the array aperture" which is defined 
as the ratio - V^/l of the collected radiant power P £ to the power 
density I in the array plane. 


Formulas for Required Nuaber of E lem e nts 

From the signal to noise considerations to be given later in the text (see 
110) we find that, for a 50A diameter transmitter antenna dish, the effec- 
tive area of the receiving antenna must satisfy the equation 


—5 D 

A eff > 2 * 10 FT “P 5 V D > W 

t T 

where P^ * transmitter power, ■ target speed transverse to line of 
sight, R x » target-array distance, * bistatic cross section of the 
target, D * diameter of the array, and A - radiation wavelength (all quan- 
tities expressed in standard MKS units). Since the effective aperture can- 
not be larger than the physical aperture, A^^ is subject to the constraint 


A < 
eff 


*D 

4 


A(D) 


(9) 


Expressions (8) and (9) require that the aperture si 2 e D be greater than 
the minimum value given by 


8 x 10~ 5 1 Vi 

n aa ■ a . 7 • > 

min w Pj. A 0^, 


( 10 ) 


The permissible ranges for ^ are then given by the shaded region in 
Figure 2 in which curve 1 is A m (r») of (8) and curve 2 is A(D) of (9). 
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FIGURE 2. ALLOWED VALUES OF A .. AS FUNCTION OF ARRAY SIZE 

ext 

2 

Bote that A^^ cannot be smaller than A m<w » v D a4n /4. Using (8) and 
(10), we obtain the following useful expression 


P 


A (D) 
A(D) 



( 11 ) 


which gives the ratio of potential thinning as a function of the diameter D. 

For highly-thinned arrays, i.e. , A^^/A «• 1, electromagnetic coupl- 

ing between the elements can be neglected. Consequently the effective 
array area will be the sum of the effective areas, A g ^, of the individ- 
ual elements. Hence, the required number of array elements (from (8) and 


(10)) is 



D . D 
. min . 
A . 


2 

v 


D . D 
min 


G 

el 


( 12 ) 


Dependence of the Required Number of Elements on the 
Radiation Wavelength. 

We note from (8) that for given angular resolution, the required 

effective receiving area is independent of the wavelength X. On the other 

2 

hand the effective area of the array elements is proportional to X jf 
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their radiation characteristic is kept constant. Hence we can conclude 

that for given angular resolution and element gain the required ainnber 

2 

of elements is inversely proportional to X . Combining (12) with (10) 
and fixing 0 » X/D, we obtain 


> 8 x 10 


-5 P T 


V 


G Tx 2 

el 


(13) 


where 



Similarly, combining (11) and (10) and fixing 6 ■ X/D, we obtain that the 
ratio of potential thinning is also proportional to 1/X^: 

8 x Iff 5 P T° 

V s <«> 

Expressions (13) and (14) indicate that increasing the radiation wave- 
length brings about two positive effects: (1) the required number of 

elements becomes smaller, and (2) the coupling between the "lenient s de- 
creases. Longer wavelengths, however, require larger transmitting dishes, 
a tradeoff yet to be evaluated. 


Some Numerical Results. 

In the following we will make use of the developed relations to find 
out what the requirements are concerning the aperture size and the number 
of elements for two different targets: (1) a low orbit satellite and 

(2) an aircraft approaching or leaving an airport. In both cases a prac- 
tical solution is sought which provides an angular resolution on the 
order of several milliradians and requires on the order of a thousand 


antenna elements. 
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EXAMPLE 1 

(a) Assuming the parameters of a low orbit satellite to be: 

5 

“ 400 Km * 4 x 10 meters 

3 

V_ ■ 8 x 10 meters/sec. 

T 

Oy ™ 10 square meters 

we find 


+8 -1 
- 3.2 x 10 sec 


which with X ■ 0.2 meters and P_ ■ 100 watts (assumed the limit to the 

t — — 

available power in the Illuminating satellite) yields (8) ■ 410 meters. 

This ' .ue implies that the angular resolution must be better than 
-3 

^max ™ 0. 49 x 10 radians and that * - integration time cannot be longer 

-3 

than t ■ 24.5 x 10 seconds. From (12) above, the number of elements 
max x2 

required (since for a filled array ■ — ) is given by 


vD' 


min 


7 ^ • -T7 “ 13.2 x 10^ elements 
4 X 


Nftin is too large for this application to be considered. 

(b) If instead of X » 0.2 meters we assume X « 1 meter, and leave 

the transmitter gain unchanged, the following results will be obtained: 

-3 

D . “82 meters, N . » 21,000 elements, 6 ■ 12 x 10 radians, and 

min min max 

T max * 0*61 seconds. The situation is considerably improved since N m ^ n 
is much smaller but is still far from being practical; for besides the 
large number of elements required, the resolution is not very good. Thus 
we conclude that imaging of low orbit satellites requires an impractically 
large number of array elements unless the illuminating power is many times 
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larger than 100 watts. 


E XAMPLE 2 

(a) For an aircraft target with parameters 

5 

• 100 Km ■ 10 meters 
Vy - 100 me ter s/sec. 

0y • 10 square meters 

we find 

Py “ 10 6 sec * 

Again assuming X ■ 0.2 meters and P^ ® 100 watts wc obtain from (io) 

®min * ^*275 meters. Suppose that angular resolution of 0 - 5 x 10 

radians Is required and that array elements with effective area A , * X 2 /2 

eJL 

(or *5.8) will be used. Then from D ■ 40 meters and (12) we 

obtain that the number of array elements will have to satisfy N > 2250 
elements. The ratio of the potential thinning is 



0.03125 


(b) If we increase the wavelength by a factor of two and require the same 
angular resolution, the following results will be obtained: D m<n ■ 0.6375 

meters, D * 80 meters, N > 563 elements, and p - 0.00797. 

On the basis of the above results, we conclude that imaging of air- 
craft can be considered practical. The required number of elements even 
for relatively small illuminating power (100 watts) is of the order of 


For instance, from (13) it follows directly that if we want to employ 
an array 400A in diameter with 1000 elements each with 7 dB gain for 
imaging of low -orbit satellites (where P T = 3.2 x 10 8 sec -1 ), then at 
X « 0.2 meters 138 Kw is required, whereas at X « 1 meter, 6 Kw is 
required. 
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1000 elements which ie reasonable. Thus It appears that microwave holo- 
graphic imag ing can find useful application In aircraft surveillance and 
traffic control systems . 

4. ARRAY SIGNALS 

signal at each array element is formed as a linear combination 
©f the direct signal and the signal reflected from all tht targets that 
during the period of observation are in the Illuminated vol ume ♦ Therefore 
we can find its properties by analysing the simplest case when only one 
target is involved, as shown in Figure 3. 

The signal measured at the point R (R is a reference point within the 
array) will be given by 

*v (t) " 8 rd (t) + 8 rr (t) (15) 


where e^Ct) is the signal coming directly from the satellite S and (t) 


la the reflected signal from the target T. 




FIGURE 3. SATELLITE-TARGET- ARRAY GEOMETRY 


ie 

The coefficients of that linear combination are determined by the an- 
gular position of the targets and the form of the radiation characteris- 
tic of the antenna. 


REPRODUCIBILITY OF J'U 
ORIGINAL PAGE IS P0< -L 



-18 


We ate Interested In the form of s r (t) when a CW signal is radiated 
from the satellite. Assuming for the radiated signal 


s (t) » CcObW t 
s o 


(16) 


the expression for e , (t) is easily found since d is constant 

ird si* 


®rd (£) “ C d cos [ w o (t ' “IT 0 ] " C d C08( V + V 


(17) 


where 


♦ 


o 



(18) 


and c is the speed of light. 

In order to find the expression of the signal s rr (t) we notice first 
.that the signal at the target has the form 


s (t) ■ C cos 

C t 



d ts Ct) -l 

“T^J 


(19) 


where is some constant. Since what is happening at the target at the 

d tr (t) 

moment t will be happening at the point R in t + we can write 


d tr (t > 

8 rr (t + > " C r cos 


[■ 




>] 


( 20 ) 


It is obvious from (20) that the signal 8 rr ( fc ) 1® a cosine function with 
time dependent phase, i.e., 
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8 rr (t) - C r cosCu o t + <Kt)) 


Substitution of (21) in (20) gives the equation 


T d tr (t) d tr (£) 

COS 0) (t + ) + ♦(t + — ) 

IOC c 


ij - cos |j 


, t d ts (t \ 

U-(t ) 

o c 


which is satisfied for 


d tr (t) d M (t) d tr (t) 

♦<t + — ) " -» 0 (-^ + — ) ♦ 


where n is integer, or for 


a i 


d £& + tr r 

; (t + . -w 0 - £ — g - 5 - 


where the dots above and denote derivatives of these functions 

T * ts tr 

with respect to t. Assuming that the change of $(t) in the very short 
time interval 


d. (x) - d, 
i tr _ _ tro i 


is negligible, where x is the observation period and d trQ - d tj .(t o 0), the 


left side of (2*',) could be replaced by 


♦ (t + 
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1* 6* | 


}<t + 

c 


-0) 


a a 

ts + 
c c 


a 


i + 


tr 


(26) 


or 


$)t) A -w 


a t8 <t - %°-> i tr <t - %s» 

~ 1 + — — 


1 + 


4 tr<‘ - 


(27) 


tr 

Since — — « 1, for any t, we approximate 


< j ) 


o,% d frp % . ; ^ d tro VN „ d tr^ “ c* * 


Kt) - - ^(d ts (t - -v*) + d^ (t - -sa )} (i - 


tr 


(28) 


For targets whose accelerations with respect to the transmitter and 
the receiver are not excessive, as is the case with aircraft and the satel- 
lites, the argument t - ^tro^ c ex P ress ^- on (28) can be replaced simply 
by t, i.e.. (28) can be approximated by 


0) 


♦ 


<d„ + K,1 U ' 


tS tiT 


(29) 


If the radiated frequency is In the L-band, then for aircraft targets with 
speeds not exceeding several hundred meters per second further simp 1 ■? f ica- 
tions is possible 


6 - - — (5 + d ) 

y C ts tr' 


( 30 ) 
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The value of the neglected term will b_ 6f the order of only 0.01 Hz. (In 
the case of satellites, application of expression (30) will lead to errors 
of the order of 10 Hz). From (30) it follows that 

♦ ■ - T d U< C > - f d tr + ♦<, < ll > 

where is some constant. The large distance between the geostationary 
satellite and aircraft allows one to consider the relative speed 


v m d 
ta ts 


as being time independent during the observation period. Therefore (31) 
can be rewritten as 


♦ - V - T d tr (t > + *o 


where 


&» ■>„ 

— v ; — — m — 
C ts* X C 


Thus, in the case of the air traffic control application, the expression 


(21) assumes a relatively simple form 


S rr (t) “ c r C0S K“ 0 + w i) c " ~l~ d tr (t) + $ 0 J 


The signal received at the ith antenna element at time t due to the di- 
rect signal from the satellite and a single reflecting target (aircraft) 
will be given by the following expression 


s 1 (t) « Acos (u> o t + <^ oi ) +B 1 cos[(w o + i^) t - y d fci (t) + ] 


Acos(u> o t + , oi ) + B^cosKi^ + «j)t - -r~| R(t) j + 1 + $ Q ] 


( 36 ) 



- 22 - 


wbere R(t) is the position of the target with respect to some reference 

element of the array and $^[R(t) ] is the phase at the ith array element 

» 

with respect to the phase at the reference element. 

5. SIGNAL PROCESSOR CONSIDERATIONS 

In this section we shall develop the main features of the spatial and 
temporal characteristics of the signal processor required for imaging of 
targets. 

5.1 Spatial Processing 

From (36) above we can start with the following expression for the 
signal received at the ith antenna element at time t due to the direct sig- 
nal from the satellite and a single reflecting target (multiple target 
considerations will be a simple extension of the arguments presented here) : 

S ± (t) » Acos(ii> o t+$ oi ) + BjCos [ (u> o +a>j) t - -y^R(t) | + ♦ li (R(t)]+4 o ]+n jL (t) 

(37) 

Here is the frequency of the carrier r sinusoid; A and are the ampli- 
tudes of the direct and reflected signals assumed equal at all array ele- 
ments; and are their phase shifts with respect to a real or virtual 

reference; n(t) is the noise present due to thermal and shot effects* 

For processing purposes it is sufficient to consider the complex en- 
velope function as follows: 

j*oi j(«.t-2i|R(t) !+♦„+♦.. [R(t)J) (38) 

S i (t) - Ae + B lC 1 X 0 11 + n ± (t) 

Note that the operations to follow can be theoretically performed equally 

well on the raw rf signals or on the equivalent if signals that can be 

envisioned at each array element. 
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Two processing schemes shall be considered here: coherent detection 

and square-law detection. 

5.1.1 Coherent Detection 

Coherent detection requires the establishment of a known phase refer- 
ence sinusoid synchronized at all array elements either at the carrier fre- 
quency or at a conv nient LO frequency. The spatial processing signal 

v (t) is formed, given by: 

• c 


v.(t,n) s E a^Ct)* -j0 i (n) (39) 

Here the symbol n is used to represent a desired "beam pointing" direction; 
o^and 6 i ( n ) represent amplitude and phase weights to be applied for beam 
formation in the desired direction and with the desired sidelobe structure. 

It can be expected that the will be approximately unity, not vary- 
ing very much with pointing direction. For optimum reception of the reflected 
signal, we point the beam in its direction, i.e. 


W * *iii R ' e >] 


(40) 


The resultant processor performance based on p 


elements can be tabulated as follows: (assuming 14 


randomly located antenna 
is input noise power) 
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DIRECT REFLECTED NOISE 

SIGNAL SIGNAL SIGNAL 


2 

Input Power A ■»F p . 

(single element) 2 



So m P, 


Nin 


Output Power a* 

(n g elements) ~2 ~ 2 ~ 


n Ho 
e 


( 43 ) 


SABLE 1. COHERENT SPATIAL PROCESSOR PERFORMANCE 


It is seen that the refleeted signal is enhanced with respect to both 
the noise and the direct signal by the factor n^. Now we can expect the 
following input inequalities to be valid (7) 




Nin 


« P 


Din 


(42) 


As a result of the spatial processing the P„/P„ situation should be ita- 
proved considerably* but the P^/P^ may still be smaller than desired. Two 
ways of enhancing this ratio would include 

a. Adaptive nulling: adjustment of ct^ and 9^(n) to produce a null 

in the resultant sensitivity pattern in the direction of the 
satellite. 

b. Doppler filtering or MTI: making use of the frequency separa- 

tion Abij' between the Direct and the Ref lee’ ' signal components 
of v c (t,n). 

5.1.2 Sauare-law Detection 

This is a direct analog of the type of processing implicit in the op- 
tical hologram. It has a great virtue of not requiring ah explicit rf 
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reference signal replicated at all the array elements. First the output 
of each array element (after rf and possibly i.f. amplification and filter- 
ing) is passed through a square-law detector, giving 

V SU lt) ■ lf 1 <t)| J -A J +2AB 1 co,{» 1 t-^|R(t)|4+ 0 ^ 11 (|EW|)-» ol } + 2*aJ 1 (t) 


+ negligible BjXB^, B^xn and nxn terms (43) 
2 

Note that the A term is a constant nv'.'J be reracr- ' essentially by 
a simple high-pass filter. Now the spatial processing can be carried out 
by forming the cos? lex signal: 


Y SL (t ’ nl) 


l 

i 


<W C) • 


. 2 , 1 - J ^ 

A )a ± e 


(44) 


from (43) the reflected signal can be enhanced, by using a 
<£(n) equal to: 


1 

i 


2 1 


± (4 u (|R(t)|) - 4 oi ) 


(45) 


The reflected signal component of the output will then be: 


V 

.sc 




j[+( U] t~|R(t) N 0 )J 

» AB,n e 
1 e 



j[*(« t- 2 ^|R(t)I+4 o -2e*(n R )J 
e 


(46) 


For a random array the second term will be negligible. The output 
noise component will now be of the form: 
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r 1 

Tsi. (t 'V» ’ 24 [ “lo < t)e 


the resultant S/H will now be: 


( 67 ) 


*Rout m 
^fiout 


A 1 

1 e 

4A 2 q M 
e o 


K 



( 48 ) 


Comparing this result with the performance of the coherent processor (see 
Table 1), it is seen that the square-law detector has 3 decibels less S/K 
in the outputs than the coherent processor. 


5. 2 Temporal Processing 

For the envisioned illuminating power, the SNR of the resolution-cell- 
outputs (outputs corresponding on a one to one basis to the resolution cells 
and coming right after the spatial processing block) will be too small to 
allow direct detection. Further enhancement is necessary and can be done 
either by using a bank of doppler filters for each output, or by taking the 
Fourier transform of each output. Both methods lead basically to the same 
result - enhancement of the SNR by a factor of 2 WTj (in the best case) , where 
2W is the system's IF-bandwidth and l/x^. is the bandwidth of the doppler 
filters. Only the Fourier transform method will be considered here. 
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Let v(t) be a baseband signal representing the output of a resolution 
cell; nultiplied by the weighting function e and let it be inte- 

grated from tine t - 0 to tine t ■ Tj as follows giving 
T 1 

V(f) - / v(t) e ~ j2irft dt (49) 

o 

Since in general v(t) is a superposition of the useful signal (target echo), 

r (t) * /P expj(2xf _t + $), and the noise, z(t) the Integral V(f ) will 
sp a 

be a random variable with mean determined by the useful signal component 
(if such does not exist the mean is equal to zero) and variance determined 
by the noise component. Thus 

^ T / j(2*f t+*) 

« v (f) - EV(f) - / /T~ e e J dt 

a p 


-j2ir (f-f d ) (a+j i sin[x(f-f d )T] 
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where (a, a+t) (see Figure 4) is that portion of the interval (0, tj) 
within which the useful signal r(t) is present) and 

oj(f) - E[(V(f) -m v (f))(V(f) ~m v (f))*] 

- E / / z(t)a*(u)e” J2wf(t - u) dtdu (51) 

o o 

• N zo T i rect( lw ) < 52 > 

where W is the low pass bandwidth of the system and N Is the power spec* 

zo 

tral density of the noise z(t) . 


1 . 


.. _i 


r 1 

0 a 


i — 

t 


FIGURE 4. TIME INTERVAL RELATIONS 


In obtaining (52) we made the assumption 


» 


1 

W 


For the signal to noise ratio defined as 

wi 2 



(50) and (52) yield 


snr t 


P T 2 


32 

N t 


zo 


I 


(53) 


(54) 


(55) 


The result (55) shows that the SNR will be maximized - the maximal 
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value being equal to 


Slffi nttx' - 5 E- 

20 


(5 '>) 


- if the integration time x^ is chosen to be equal to x i.e., if we start 
integrating exactly at the moment of appearance and end the integration 
at the moment of disappearance of the useful signal. Shorter integration 
time leads to seller SNR since we do not make use of the whole available 
energy in the target echo. Longer integration time, on the other hand 
leads again to smaller SNR since in the excess time x^ - x we integrate 
only noise and not useful energy. 


5.3 Block Diagram of the System 

As was mentioned in 5.1, we have considered two processing schemes; 
coherent detection and square law detection. We found that both are appli- 
cable to the system. Coherent detection has the advantage of 3 dB greater 
SNR, however the latter has to be obtained by establishment of a known phase 
reference at all array elements - a requirement which is not easy to ac- 
complish in a relatively large array of hundreds or thousands of elements. 
Square law detection does not require an explicit reference signal repli- 
cated at all array elements. This is a great virtue and was the reason 
for giving priority to square law detection in almost all power level 
calculations that have been done throughout the text. 

We have come to the conclusion that the most adequate form of pro- 
cessing for the system is the digital signal processing and that, in order 
to meet the requirements of practical application, it has to be done in 
parallel, i.e., the res "ution cells of a given sector or of the whole 
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volume of interest should he interrogated simultaneously. 

Except for the differences in the array element modules, the system 
block diagram is basically identical for both processing schemes. A 
simplified representation of it is given in Figure 5. 



FIGURE 5. SIMPLIFIED BLOCK DIAGRAM OF THE SYSTEM 

The array outputs are sampled at a rate of f ■ 2W samples per second, 

8 

where W is the bandwidth of the system in baseband. If one set of samples 
is visualized as an N-dimensional vector _h, where N is the number of array 
elements, then, the set of subsequent vectors sampled during the Integra-; 
tion time ij will constitute an NxM-dimensional matrix of the form 

M « 2WTj 


I h^ h ^ lit » 


(57) 
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Thls matrix is the processor Input data. * The output of the processor can 
be represented by an N.. xM 1 -dimensional matrix l.e.. 



( 58 ) 


where la the number of the Interrogated resolution cells and is the 
number of Fourier transform coefficients. Each of the rows is a discrete 
Fourier transform of the signal associated with the corresponding resolu- 
tion cell. 

The data processor will have to digest and interpret the output data 
and prepare for display only condensed information. 

Figure 6a shows the block diagram of the array element module for 
square law detection. 


I PROCESSOR 
* f c s 2Wf 



FIGURE 6a. BLOCK DIAGRAM OF ARRAY ELEMENT FOR SQUARE LAW DETECTION 

The block diagram of the array element module for the case of coherent de- 
tection, together with the block diagram of the reference phase generator 
Is shown in Figure 6b. 




FIGURE 6b. BLOCK DIAGRAM OF ARRAY ELEMENT FOR COHERENT DETECTION 
AND REFERENCE PHASE GENERATOR 

In both cases, the direct signal component Is filtered out In baseband; 
for the coherent detection case the filtering could have been done In the 
IF-band, as well. 

5.4 Nulling of the Illuminating Source 

The system will have to detect reflected target signals which are 
as weak as 100 dB below the level of the direct ilgnal. The transmitter- 
noise power measured in a 1 Hz bandwidth (usually 80 to 130 dB/Hz below 
the level of the carrier, depending on the kind of transmitter and the 
distance in Hz from the carrier) is comparable with the target reflec- 
tions. Further, the significant components in the transmitter-noise 
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sidebands extend beyond any conceivable doppler frequencies [ 5 ], This 
means that, unless suitable measures are undertaken, we may start "de- 
tecting" false targets in various directions and with different "doppler" 
frequencies - a result of the penetration of transmitter-noise componentJ 
through those sidelobes which happen to be in the direction of the 
transmitter. The problem of transmitter noise can be avoided if we atten- 
uate the direct signal by suppressing the array response in the direction 
of the satellite. A simple way to do this, provided the position of the 
illuminator is known, is to use array elements with sensitivity charac- 
teristics having nulls in the direction of the illuminator. However it is 
obvious that such a technique works only in the case of synchronous defec- 
tion, it doesn't work if square law detection is employed since, in the 
latter case, the presence of the direct signal is pertinent.. Nevertheless, 
the idea of using arr-'v elements with nulls in the direction of the illumi- 
nator will help us to develop a nulling circuit which could be applied when 
square law detection is used. 

Square law detection requires not only presence of the direct signal 
at the input of the square law detectors, but the direct signal must be 
at least several times stronger than the corresponding thermal noise. It is 
clear from the expression (59) describing the output of the square law de- 
tector (hero A, JJ and Z are complex amplitudes of the direct signal, a re- 
flected signal and the narrow band thermal noise) that if the direct signal 

V« (|A | 2 + |fl | 2 + |z| 2 ) + (AB* + A*B) +(AZ* + A*Z) + (|z | 2 -|z| 2 ) 

(59) 

Constant Terms Signal Terms Linear Noise Quadratic 

Term Noise Term 


[5] 


K. W. Saunders in M. I. Skolnik (ed.), 
McGraw-Hill, Inc., New York, 1970. 


"Radar Handbook , 1 Chester 1, 
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1s suppressed too much (i.e., below the level of the thermal noise) then, 
because of the quadratic noise term, the SNR will no longer be independent of 
the direct signal but will be proportional to it. In other words, suppres- 
sion of the direct signal before performing square law detection may mean 
a decrease cf the SNR, which cannot be tolerated. Therefore if the null- 
ing of the direct signal has to be done, it should be done aft^» the square 
law detectors. 

In the following, two methods of nulling - one for a system with syn- 
chronous detection and the other for a system with square law detection - 
are described. Both lead to array elements with circularly symmetrical 
radiation characteristic suitable for application in an air-traffic sur- 
veillance and control system. 

We will pursue the idea of usin^ array elements with nulls in die 
direction of the geostationary satellite. Clearly, it is desirable that 
>we should be able to adjust the pos ' . ■ the null whei: necc t *ry. Not 
only should we not expect the geosynchronous satellite to be strictly sta- 
t. onary, but we must also allow for the possibility that it can change con- 
siderably its longitudinal location, if, for different reasons, such a 
change is required. 

Consider the circuits in Figure 7. In both cases it can be shown 
th/’t the deterministic part of the output will have an amplitude propor- 
tional to: 

l s out (0) l " l G(0) l 1 F(0) I (60) 

where G(0) is the individual element sensitivity characteristic and 

| F (0 ) | * 2|sin(^l (sin0 -sinft))| 

X 


( 61 ) 
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FIGURE 7a. ARRAY ELEMENT GEOMETRY, 
SCHEMATIC FOR SYNCHRONOUS DETECTION 



where is the elevation angle of the satellite srurce. In the synchro- 
nous detection case it is necessary to set the angle 4> to the value: 

$ « sinC (62) 

In the square law detection case, the corresponding adjustment is 
done automatically. 

We must note here that the expression (60) has been obtained under 
the assumption that the sensitivity characteristics of the dipoles and 
D 2 and the input-output characteristics of the square law detectors SLD^ 
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and SU >2 were identical. Since in practice the characteristics of the ele- 
ments differ from element to element , the^e will be leaking of energy for 
0 , i.e. , the nulling will not be as good as indicated by (60) and (61) . 
Still a 10 dB or greater suppression of the transmitter noise can be ex- 
pected, if the elements in the pairs constituting one array elenent are 
carefully chosen to have as close characteristics as possible. 

6. ANGULAR RESOLUTION AND FOCAL DISTANCES OF A CIRCULAR ARRAY 

It can be shown that in the vicinity of the focusing point the radia- 
tion pattern of a random array approaches the radiation pattern of the 
corresponding contfnuous aperture, the variance being equal to 0.5/N at 
the half power points (N is number of elements) and going to zero at the 
focusing point. For ease of analysis the continuous circular aperture is 
taken as the model for the random circular array. It will be shown that 
the near field azimuthal resolution of 100-wavelength or bigger arrays - 
for distances as small as R * 40 p q , where p q i„ the radius of the array 
aperture - is approximately given by the same expression that applies for 
far field conditions. Formulas for the angular resolution in elevation 
for far field conditions will be presented. The array focal distances 
for broadside and endfire pointing will be determined. Although most of 
the time we will find it appropriate to use the term "continuous aperture" 
we preserve the freedom to speak about "array aperture" or simply "array" 
when convenient. 

6.1 Angular Resolution in Azimuth of a Circular Array 

Consider a circular array with radius focused at the point F and 
illuminated by a point source T. Let, as shown in Figure 8, the positions 
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of F and T differ only it :imuth, i.e. , Pp * P T * R Q » °p " °x " ® 0 * 

“ 0, « Af. We are interested in the output of the array as a func- 

tion of the aaiauthal difference &$. 



FIGURE 8. GEOMETRY FOR AZIMUTH RESOLUTION CALCULATIONS 


The CFA (complex field amplitude) at arbitrary point M(p, $) in the 

array plane will be given by s * Ae TM , where r^, is the distance be- 

TM 

tween the point M and the source. The complex weight associated with M is 
j kr 

w * e J FM, where r^ is the distance between the focal point F and the 
point M. Since the array output is weighted integration over the array 



pdpd$ . 


(63) 
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Expressing r ra and in terms of the coordinates of P, T and M one obtains 

2 2 ?< 
r IM “ t R Q + P - 2 R o p cos« 0 cos (4 - A<j>)] 

r FM “ {R o + p2 " 2 V coso o (64) 

In the discussion to follow we will assume 


» 50> 
o 


(65) 


Also, we will confine our interest to only small values of A$, 8ay (55) 


t 0 


( 66 ) 


and to ranges R q such that 


R > 40p 
o o 


(67) 


Under the conditions (66) and (67) the difference can be 

fairly good approximated by the first two terms of its Taylor series ex- 
pansion in powers of p, i.e., 


r TM “ r FM “ cosa o (cos<J> - cos (<J> - A <fc) )-^y 


cos 2 a 


R — (cos 2 $ - cos 2 (4> - A4 >))|t 


( 68 ) 


Making use of the approximations sinA$ * A4>, cosA l • 1 and ignoring (A$>)' 
we obtain 


r TM " r FM A " (cosa o sin4,A^)p 


- (cos 2 a sin2$ Ai}>)— 
2R 


(69) 


x ur ItiP 
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Sub8titutlon of (69) In (63) yields ^ 

P o 2 k J (kpcosa A*)sinf J (k 2R coa\ti)mia2* 
v(A$) -a/ / e 0 e ° 


o o 


pdpd$ 


But, since (66), (67) 


cos 2 

o 


a Ad 
o T 


2k 

X 


2R 



one can approximate 


(70) 


(71) 


j (k|^co8 2 a o A6)sin2() 

e ° A 1 + j0t^-co8 2 a o A^)sin2^ (72) 


Thus 

p o 2k j (kpcosa A$)sin$ 
v(A+) -4 J e ° pdpd$ 

o o 


j R 

'2R P o 2k j (kp cos a A$)sin$ 

+ 3 °cos 2 a A^) / / e 0 Bin2*p 3 dpd* (73) 

0 oo 

The second integral is however equal to zero, as can be easily seen if we 
integrate first with respect to and we have 
p o 2k j (kpcosa A$)sin$ 

v(A$) ±4 / e 0 pdpd* (74) 

o o 

Hence, under the conditions (66) and (67) the contribution of the quadratic 
term of the difference r^ - r RI (68) to the integral is negligible. Thus 
the array pattern for targets near the focal point will be the same as for the 
array when focused at 00 . 

We solve 


v(A$) - 2 k a/ J (kpcosa A<p)pdp 
o p 


kp 2 A- 


2J ^ (kp o rosa o A<|)) 

^ 1 ■ 1 ■ "in 

kp cosa Ad> 
o o 


( 75 ) 
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where J (x) and J.(x) are Bessel functions of zero and the first order 
o 1 

respectively. 

Figure 9 shows the plot of the function 2J^(x)/x. As could he seen 
the values of x at which 2J^(x)/x * /2?2 and 2J^(x)/x ■ 0 are denoted by 

*1 V 



FIGURE 9. PLOT OF 23^%) /x 

From 

x, “ 3.8317 ■ kp cosa 64 

we find that the azimuthal resolution <5 determined according to the 
Rayleigh criterion, is given by 

U - - 1 -^- (76) 

• V R p cosa Dcosa 

0 0 0 

where D 9 2p is the diameter of the aperture. 

0 

The 3 dB bearawidth 6 d> (measured in azimuth) follows from 

» 1.6 = kp o cosa Q • 
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Thus we have 


«♦ 


0.51 A 

pCO80t 
O 0 


1,02a 

Dco3a. 


( 77 ) 


Froa the discussion above it is apparent that the expressions (76) 

and (77) can be used for distances even smaller than 4 Op . (40 is not a 

o 

aagic number) if of interest are only approximate values. 

2 

Table 2 compares the near field — far field transition range Rj“8p o /A 


and the distances R 


40p for four arrays, 
o 


A » 0.2 (ffl) 

Rj^Cm) 

R 0 (m) 

P o - 50A 

4000 

400 

P D - 100A 

16000 

800 

p 0 - 200A 

64000 

1600 

P 0 ‘« 500A 

400000 

4000 


TABLE 2. DISTANCES REQUIRED TO APPLY SIMPLIFIED 
FORMULAS FOR FOCUSED BEAMS. 

6.2 ANGULAR RESOLUTION IN ELEVATION OF A CIRCULAR ARRAY 

Suppose, again that our array is focused at the point F, but now il- 
luminated by a point source T which position differs from F only in the 
elevation coordinate by a small angle Act, as shown in Figure 10. 


The output signal will be given by 

,, . f -l k < r Df r rn ) . , 4 

v(Aa) * J J e pdpd* 

o o 


(78) 


where 


r TM " < R o + p2 ** 2 K o P co s(a 0 + Aa)cos<J>] is 

r~ . " [R 2 + p 2 - 2 R pcosa cos*] J 
FM o 0 0 


( 79 ) 
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' are the distances of T and F from an arbitrary point M in the aperture plane. 


FIGURE 10. GEOMETRY FOR ELEVATION RESOLUTION CALCULATIONS 



The integral, above, cannot be solved in a closed form and an approxima- 
tion of the form in (72) doesn't help. Therefore in the following we will 
give some results which were obtained by ignoring the quadratic phase 
factor in (80) i.e., by assuming far field conditions. 



-43- 


The angular resolution in elevation, determined according to the 

Rayleigh criterion is solution of the equation 

2 , (%>* 

t p 0 ( s v lo “o + — r~ coso o ) ■ »•«« (81) 

Hie 3 dB beamwidth (measured in elevation) can be determined from 

2» 1 C^a) 2 

T p 0 I( 2 6a) Sinot o + ~Z cosa o 3 " 1 * 6 ( 82 ) 


For the special cases characterized by elevation angles a Q » 90 and 
0°, we easily obtain: 


1.22A 

9 ■' 

0 

1.02 A 
D 


6a * 2 


1.56 /j^~ 

2 . 02 /| 


RFi'Kol^CTBlUTY OF THE 
ORIGIN Al. PAGE IS POOR 
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6.3 FOCAL DISTANCES 

Consider a continuous circular aperture with radius p Q focused at 
the point F and illuminated by a point source T which position differs 
from that of F only in the range coordinate, i.e., * R T p p =• Ry, 

“ Op ■ a Q , » 0, as shown in Figure 11. 



We are Interested in the form of the expression 

0 2ff — j k (^rKf-^ru) 

vCR^lRy) ■ / / e ^ pdpdifi 

* * ft A 


where 


r TM " <4 + p2 " 2R T Pcosa o cos*)»i 
r FM " (R F + p2 “ 2 V C0S V 0S ^ ?S 


(85) 


( 86 ) 
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are the distances of T and F from an arbitrary point M(p,4») in the aper- 
ture. Note that the array sensitivity along the direction of focusing is 
given by 

P 

sagRp - x ivogRp)! 2 (87) 

T *T 

where: P g * array power output and - radiated power from the distance 


Expanding the difference r ^ I “ r FM In a Taylor series in powers of 
p and ignoring all the terms beyond the quadratic term, one obtains: 

r M- r ra" II T- R F + ^-^> (1 - co„ 2 « 0 • co. 2 «f 2 (88) 

Letting 



The last integral can be solved in closed fore for two special cases 
of the elevation angle: (a) for o q ■ 90° and (b) for ■ 0°. The re- 

sults for these two extreme situations give a fairly good insight of what 
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ia happening for other values of a^. 

In the case (a) (broadside focusing) the Integral (90) assumes the 


v(a) « 2 v £ e pdp 


and we obtain 


ap 2 ap 2 

. 0 , , 0 \ 

"j 2 sin ( 2 


or substituting for "a" from ( 89 ) 


2 r£rJ: I'** 2 ! 


1 1 sin 2 [-t-C y — p- 

|v(a)| 2 - |v£- f-)| 2 « 4 h 

“t F 1 l x ,, 




In the case (b) (end-fire focusing) the integral (SO) assumes the form 


P A an Z 


v(a) 2 ir/ e J C-5 2 -) 


f / *.<«> * 


j rt ap- ap- 

•»o • < J o (_ r ) + J J i ( i a)> 


where the last integral was solved by using the relations 


2 4- .. Jz v+1 

/ A\(tMt . -s-^rr 


2 v+1 (J v (z) ‘ J J v+l (2 » 
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for v ■ 0 [6]. Combining (93) with (89)’ and squaring both sides we obtain 
the expression 



J o [ 4 C P^ " R p )p c J * + J l l 4 ( R^ " R^pjl 


(94) 


2 2 2 

The curves of the functions f, (x) = (sin/x) and f-(x) ” J fc (x) + J, (x) 

1 2 o l 

(both f^(x) and f^(x) are even) are depicted in Figure 12. 

By aid of tables we find that the half power value of x for the first 



a. Broadside b. End fire 

FIGURE 12. FOCUSING ARRAY SENSITIVITY FACTORS FOR 
BROADSIDE AND ENDFIRE POINTING. 


function (case a) is x . ■ 1 1.395 whereas the half power value of x for 

the second function (case b) is x _ ■ i 1.69. Employing the equation 

0x5 


— p2 / A _ A \ 
4 P o Rp ; 


= X oB } 


[6] "‘Inndbook of M ithimatical Functions with Formulas Graphs and Mattie 
critical Tables" , National Bureau of Stand'-du, Applied Mathematics 
Series 55, p. 433, June 1964. 
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• £ 

we easily obtain that the hyperfocal distances 
** „ 

distances are given by 


and the near-focal 

2 

case a. =■ D /3.55A 

case b • 

^Hb 

- D 2 /4.31A 



Sib 

» D 2 /3.6?A (96) 


where D “ 2p^ is the diameter of the array- Note that the focal distances 
are about 20Z smaller in case b (endfire focusing) tuan in case a (broad- 
side focusing). (Examples With A *° 0.2m and D * 400A “ 80m, R, “9 km, 

da 

Rgk a 7.43 km and R^ a = 4.5 km, « 3.715 km), for 0 < a Q < 90° the 
focal distances lie between the values g .ven in (96) . 

It is not difficult to see from (95) that if the array is focusad at 
the distance Rjj/k, where k is an integer, then the distances at which - 
ignoring the inverse square spreading loss - the array sensitivity decreases 
3 dB are given by S H /(k-l) and R^/(k+l). It follows, therefore, that if 
a 3dB loss is tolerable and if infinity is chosen as one of the focusing 
distances, then the other focusing distances can be determined from the 
relation 

Rp n “ Rjj/2n (n = 1,2, ... ) (97) 

We will call the distances defined by the las', expression - including in- 
finity - as ”3 dB focusing distances”. 


The hyperfocal distance is defined as that focal distance at which 
the array is in focus outto infinity, i.e., the power density at R , 
ignoring inverse square spreading loss, is 3 dB above the valea at h 
infinity. 

The near focal distance R^, is defined as the range on the near side of 
the hyper. focal distance all which the power density, when the array is 
focused at the hyperfocal distance and the Inverse square spreading law 
is neglected, is 3 dB less than at the hyperfocal distance. 
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• 7. DEPENDENCE OF THE POSTINTEGRATION SNR ON THE RANGE AND THE SPEED 
OF THE TARGET WITH THE INTEGRATION TIKE AS A PARAMETER 

7 . 1 Introduction 

The targets that are to be detected move with different speeds, at 
different ranges and in all possible directions. As a -esult the echoes 
appearing at the beam outputs will have different durations and different 
amplitudes. 

The pred? tection integration performed by an FFT algorithm will assoc- 
iate with these pulses, and therefore with the corresponding targets, post- 
integration signal- to-noise ratios with different magnitudes. It is con- 
venient to use this postintegration SNR as the basic measure of system 
performance. The effects of further filtering of the detected signal will 
be studied later. In other words, the post integration SNR is function of 
the target characteristics such as bistatic cross-section, range and trans- 
verse speed. Accordingly, for assumed bistatic cross-section, given dura- 
tion of the integration interval and given array and transmitter charac- 
teristics, one c>n represent the post integration SNR as a surface above 
the "R-V plane" where R is the range and V is the transverse speed of 
the target. The part of the R-V^ plane above which the associated SNR ex- 
ceeds a required threshold level can be viewed as "domain" or "re 0 jen of 
coverage" of the system. 

Our objective here is to determine the • iape of this region and its 
dependence on the duration of the integration interval. The obtained re- 
sults will answer the following questions: (1) what is the required power 

to cover a prescribed region in the R-V plane - given the integration 
time, and (2) what is the integration time for which the required power 
is minimized? 


reproducibility of the 
ORIGINAL PAGE IS POOR, 
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In the discussion to follow we will assume that (1) the boundaries 
between the beams are sharply marked and that therefore the duration of the 
target echo is determined by 



where V is the target speed and d Q is the length of the target path inside 
the beam; (2) that both the frequency and the amplitude of the target echo 
pulse appearing at a beam output, are constant. Although this means a great 
deal of idealisation the results obtained will be useful to conceptualise 
the system requirements. 

Later, (Section 7.4) by taking into account the change of the doppler 
during the integration, we will extend the discussion to the cases for which 
the assumption of constant frequency pulse doesn't hold, and thus obtain 
more general and of course more accurate results. 

7 . 2 Postintegration SNR 

In approaching the problem of the post integration SNR we will assume 
that all beam outputs are repeatedly integrated every seconds in such 
a vay that no portion of time is left outside the integration intervals. 

In other words if we start integrating at t * 0 we will end the integra- 
tion at t ■ and immediately start a new integration interval. Our at- 
tention will be concentrated not on the results associated with a particu- 
lar integration, but on the results associated with the passage of target 
* 

through a beam. Let us concentrate on a beam output in which a target 

We emphasize this fact since, as can be understood from the above, another 
approach is also possible where priority would be given to the results fol- 
lowing a single integration interval. More precisely, in this approach it 
would be required that a t leas t one of the beams that are run through by the 
same target during an integration interval would give sufficient bMR. 
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echo appears during one of the integration intervals and let for convenience 
the interval of appearance and the intervals coaing after it be denoted by 
1^, I^j and so forth, as shown in Figure 13. 



, ‘ t: ^ 

i . ^ ^ 1 

1 

S . 

• 1 • 

! „ T, ! | 

1 

» 

1 

1 

1 

I 

1 L_ 1 

1 

1 Tr 

! t s ! t e J 

1 1 1 

« I 1 

0 Tj 

t 

1 

1 

3TJ 

! t 

1 

1 


FIGURE 13. TARGET ECHO INTEGRATION INTERVALS 


The target echo, which can appear at any moment t « t fl of the interval 1^, 

disappears at the tine t = t + x in the same or in one of the next inte- 

e s o 

gration Intervals, depending on its length t^, the length of the integra- 
tion time T_ and the starting moment t . Obviously, here we cannot talk 
A s 

about postintegration SNR in the same sense as we could if the target 
echo were completely contained in the interval of integration. The trouble 
is that for the same target echo we can actually have more than one post- 
integration SNR (one for the interval 1^ and another for the interval 
in the case of Figure 13) and that, in addition, each of these SNRs can be 
considered as a random variable. (We can view the set of values that can 

be assumed by the time t g , i.e., all the points belonging to the interval 

Ij, as being the space of definition of these variables). 

To avoid this problem it seems reasonable to define a bound on the 
postintegration SNR as follows: 



(98) 

Namely If Tj. < x^/2, then at least one of the integration intervals will 

be inside the target echo interval, and according to (55) the SNR associated 

x 

with it will equal P x_/N ; on the other hand if x_ > -r- , then one of 

sp 1 o 12 

the integration intervals will overlap with at least half of the duration 

of the target echo and according to (55) the corresponding SNR will be 

equal to P (x /2) 2 /N x, . So defined, SNR implies design for the "worst 

sp o o I 

case" so that equation ^98) will be a lower bound on the SNR achieved. The 
dependence of the postintegration SNR on the integration time is illustrated 
in Figure 14. 



FIGURE 14. DEPENDENCE OF THE MINIMAL SNR 
ON THE INTEGRATION TIME 

We see that MSNR assumes its maximum value at x^. ■ t q /2, when the in- 

tegration time is exactly one half the duration of the tar -*' ho. Note 
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that the maximum value of the MSNR is only 50% of the theoretical maximum 
given by 56. 

7.3 Range and Speed Dependence of the* Post integration SNR in a 
Microwave Holographic Imaging System 

In this section we will develop expressions in a more practical form 
in terms of parameters like the size of the array aperture, the target range 
and the target speed. Expressions which describe the dependence of the MSNR 
on the target range and target speed will be obtained. It will be tacitly 
assumed that we are thinking in terms of a system where only low elevation 
angles (say up to 20 ) are of interest. Making use of the bistatic radar 
equation (53) and considering the case of square law detection one can write 

P - J- . p t „<r A e ;t . _ L (99) 

5P ~ 2. 4-^ R> 

where 

P ° transmitter power 

G t * gain of the transmitting antenna 

R ■ geostationary satellite - receiver distance 
8 

0 « Bistatic cross-section of the target 

R * target-array distance (100) 

The factor 1/2 appears because half of the power is lost in the 
"conjugate image" substitution of ( 99 ) and 


XR 
" X>v t 


( 101 ) 


in (98), where D is the size of the array aperture and V is the tangential 
(azimuthal) component of the target speed yields 



MSHR (R,V. ) - 

w 
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<^? t a-tKo^_ XR 

2M 2 N 0 Rs « z * 2P ' t i. 


| ^ t P t rAef( i 

[a(Atr)X « s i%4 F'rjV? 


V* £ 


2DT t 


ao2> 


Thus we have obtained an expression which describes the functional de- 
pendence of the M3NR on the target range and target speed for fixed in- 
tegration tine Tj. In the case of synchronous detection, the right side 
of (102) will have to be multiplied by a factor of two. 

The functional dependences of the MSHR on for fixed R ■ and on 
R for fixed » V f1 , are plotted in Figure 15 and Figure 16. A three- 
dimensional illustration of the surface (102) is given in Figure 17. 



FIGURE 15. TANGENTIAL SPEED DEPENDENCE OF THE 
MINIMAL SNR FOR FIXED RANGE 

It is important to see that all the points (R, V t > in the R - V plane, 
satisfying the inequality, 

MSNR(R,V t ) >_ MSNR q , 


(103) 



FIGURE 16. RANGE DEPENDENCE OF THE MINIMAL 
SNR FOR FIXED TANGENTIAL SPEED 



FIGURE 17. MINIMAL SNR SURFACE OVER THE RANGE 
- TANGENTIAL SPEED PLANE 
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where MSNR & is an arbitrary constant, lie in the rectangular region 

(R < R , V < V ) whose distensions R and V are solution of the set 
— o t — to o to 

of equations 


Aetf* 6 t *P t • <r ^ 

2(4if) z N 0 R S * * Zo M5NR a = R 


(104) 


V t 


2 


2DT-, 


R 


(105) 


Consequently the domain of the system, i.e.» the set of all the detectable 
ranges and tangential speeds, can be determined directly from (104) and 
(105) by simply replacing the constant HSNR^ with the lowest tolerable 
value for the MSNR. 

The following interesting conclusions can be drawn from Figures 15-17: 

(1) if we increase the range and keep the speed constant, i.e*, v t = v tr 
the MS NR will be constant until we reach the distance R * ^Dt^V^/X; after 
this point it will start decreasing inversely as the square of the distance 

(2) similarly, if we increase the tangential speed and keep the range con- 
stant i.e., R - R^, the MSNR will be constant, until we reach the speed 

V * XR^/2DXj.; after that it will start decreasing inversely as the square 
of the speed; (3) a target being at the maximum detectable distance, i.e., 
at R * R P r °duces the same MSNR irrespective of its tangential speed 
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as long as its speed is smaller than the maximal detectable speed V^. mnv » 
and conversely, a target moving with maximal tangential speed produces the 
same MSNR irrespectively of its range *as long as its range is not greater 


By changing the integration time (and keeping the other parameters 
constant) one can adjust the system to different sets of targets as shown 
in figure 18. Note that larger integration time extends the sensitivity 
of the system towards larger ranges whereas smaller integration time 
extends its domain toward larger speeds. 





FIGURE 18. INTEGRATION TIME DEPENDENCE OF THE RANGE- SPEED 
COVERAGE FOR FIXED ILLUMINATING POWER. 


The ratio of the integration times corresponding to case a and case b in 
the Figure, i.e., Xj a an d will be given by 

TlO. CL 

Txb \ ^-vnax b/ 



\ / 

v t max a. 


( 106 ) 
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From (102) we also find that the required power to detect a target at range 


E ■ E 

max 


raw ^5 ^ 

r t rmn — — 


1 ^2. 
mom 


Ae(f G t CT Tj. 


(107) 


x R 

and is Independent of Its speed If * also the re< J uired 

power to detect a target with speed V t " ^ tmax wiH be given by 


p _ 8 («*> (msnr)^ NoR tftj V t 

*tmw) “ 


V / ^ 

ma* 




(108) 


^ ^ 2.T) 

and will be Independent of its range provided R < • 

If we want the system to cover a region defined by R <_ and 

V < V (see Figure- 19), we will determine first the integration time from 
t — tmax 




X Rmox 

2SV tM , 


( 109 ) 
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»r>d for this integration time the required illuminating power will have to 
satisfy the inequality (104). 


p ^ 4 wyj)t 

1 " A e(( G, t > r (no) 

Any other choice of the integration time will require greater illuminating 
power than that given by (110). This is understandable since with integration 
time different from that given by (109) we can cover the shadowed ax«.a only 
by a larger rectangle , and therefore by greater illuminating power (see 
Figure 18). Note that for given by (109) both (107) and (108) are equal 

to the right side of (110). Smaller x_ increases the value of (107), where- 
as bigger Tj. increases the value of (108). 


7.4 Frequency deviation of the reflected signal as a factor which limits 
the length of the Integration time . So far we have been neglecting the 
effect of doppler variations upon the SNR and the length of the integration 
time. This enabled us to concentrate on those factors which are crucial 
in determining the shape of the SNR surface in the R -V t plane and the 
optimal integration time (i.e., the one that provides sufficient SNR 
with minimal illuminating power); namely the array size, the range and 
the range rate vectors. 

In thi6 section we will examine how the fact that the reflected 
signal changes its frequency during the integration period affects the 
SNR and the optimum integration time. For the sake of simplicity of the 
discussion only linear changes of frequency will be considered, i.e., the 
target echo will be described by 


*(t) = iTp 


j(Vrrfdt + 1ft*) 
e 


(in) 
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where is the main doppler shift and fX is the rate of change of the 
radian frequency. 

First we will determine the output SNR of an integrator as a function 
of the integration time if the input signal is 

v(t) 4- *(t) (112) 

where r(t) is given by (111) and z(t) is complex noise. This will lead 
to a diagram describing the degradation of the SNR a? a function of the 
rate of change of radian frequency and will give an expression for the 
"critical" integration time, i.e., the integration time which for given 
rate of change of the radian frequency maximizes the SNR# Next we will 
assume two points, one representing a stationary receiver placed at the 
origin of a coordinate system and the other representing a moving source 
which radiates a signal with constant frequency. From geometrical 
considerations we will find out what is the rate of change of the radian 
frequency for the signal detected by the receiver. Then a 

relationship expressing the critical integration time directly in terms 
of the target range and range rate vectors will be obtained. On the basis 
of these results we will be able to describe the effect of the frequency 
deviation on the SNR. Finally we will determine the optimal integration 
time and the required illuminating power. 

The analysis x^ill be based on the assumption that of interest are only 
small elevation angles and targecs moving in horizontal lines as can be the 


case if the system is used for air-traffic surveillance. 
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7.4.1> Signal to nolae ratio . Suppose the signal given by 112 and 111 
Is Integrated from — T/2 to T/2, i.e.. 


v({ It) = 


v (t) e it 




(113) 


The expectation Jtj of V ^ | ^ will be given by 


>«*(<!*) = 


Ev(f|r) 


j(2Tff a t -jeiT-ft 


(114) 


The variance 


tfCl'f) ° £ V(f \r) will be given 




= e{Cv(f M -At |T)] [v(f |r) -yu(f lr)J*J 
. f r/z r\^r. . 

-J I Z(t) Z (u)e dtda 

/V. V J 




(115) 


For the output SNR defined ns 


snr(t) - 



■(ilT) 


( 116 ) 
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(note that is the frequency associated with v.he singal r(t) right in 

the middle of the integration interval; it can be shown that /l (i | c) 
assumes its maximal value at f “ f ,) one obtains 




r‘ C/ Z 

CosAt^t + 
2 . / 



- fcPlfal 


IS 


N, 


C 2 (u) +■ 5> 2 ( u ) 


u 


. here 


and 



( 117 ) 


( 118 ) 


019 ) 


are the Fresnel integrals. 


o 
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The dependence of the SNR on the integration time is illustrated 
in Figure 20. As can be seen the output SNR is maximized for u * 1, hence 



is the "critical" Integration time. 


( 120 ) 



FIGURE 20. SNR DEPENDENCE ON THE INTEGRATION TIME. 

The diagram in Figure 20 which represents the dependence of the SNR on 
the integration time T for fixed AA can be used to obtain another 
interesting diagram showing the functional dependence of the SNR on the 
rate of the radian frequency sweep yU for fixed % . Such a diagram is 

given in Figure 21. 


REPRODUCIBILITY OF THE 

ORIGINAL PAGil IS POOR 
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FIGURE 21. SNR DEPENDENCE ON THE RATE OF CHANGE OF 

RADIAN FREQUENCY FOR FIXED INTEGRATION TIME 


The expression 


( 121 ) 

is nothing but the ratio of the SNR obtained by the integration of a 
signal with rate of the radian frequency sweep, JX , during time X 
and the SNR which would have been obtained if there was no frequency 
sweep, i.e., for /*■ - 0. Therefore it gives the degree of degradation 
of the SNR due to the variation of the frequency, and will be called 
hereafter "degradation ratio". 


F(u) = 


SNR 

P% 


Wc 
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7.4.2. The rate of change of the radian frequency of the signal received 

froa a aoving source . Our next objective will be to determine how the 

rate of change of the radian frequency of the sin gal received from a moving 

source is related to the position and velocity. For that purpose let us 

assume that ve have a stationary receiver placed at the point R g (see 

Figure 22) and a source which at t * 0 passes through the point T . 

8 



R s 


FIGURE 22. RECEIVER - MOVING SOURCE GEOMETRY 


Let, as shown, the vector position of the point T be deonted hy R , the 

8 O 

the source velocity vector (we assume that it is independent of time) by 
V^ and the source position at arbitrary moment t by R. By definition 


dcP d-co ^ 

/*■ ~ At &t 


( 121 ) 
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where w ^ is the doppler frequency shift of the signal, i.e.. 




MiJSi 

X dt 


Hence, for t * 0 





(122) 


(123) 


Since 


fal =Ir 0 +v, 


W = (| 


— 1 2 


Roi + 2R 




|V 0 |VJ 


% 


(124) 


(123) yields 



I ft/ KF- (<Vv/ 

m 3 


(125) 


For targets moving along horizontal linec expression 125 assumes the form 


/* = - 


2Tt 

X 


V<? 

Ro 


(l — COS>V COS , 2 Qj 


( 126 ) 
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uitere: V^m |V q | , R q ■ |R q |, ■ elevation angle, 9 - angle between the 
projections V op and R Qp of V q and R q on the horizontal plane as sketched 
in Figure 23 


\ 

\ 



\ 


FIGURE 23. PROJECTION OF THE MOVING SOURCE - 

RECEIVER GEOMETRY ON THE HORIZONTAL PLANE 

In the case when only small elevation angles are of interest, cos ot can 
be approximated by unity and (126) simplifies to 

i* = _ ^ «n*S 

^0 ^ Ro (127) 

where 

V tc = V o I 6 


( 128 ) 
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is the tangential speed of the target. Applying 120 we find that the 
critical integration time for this case is 

_ \|2XR 0 

c " v to < 129 > 

Note that the relation 129 associates with each point of the R - V plane 
a critical integration time, i.e., it tells us how long to integrate a 
signal coming from a source at range R with tangential speed in order 
to obtain maximal SNR. 

7.4*3* Representation of the degradation ratio in the range speed plane 
of the system * As can be understood from the above, any echo signal due 
to a target at range R and with tangential speed V^, will be chart cterized 
with the rate of change of radian frequency given by (127) . 

On the other hand, as will be shown soon, the integration time of the 
useful energy, X » which accounts for the minimal SNR is also a function 
of R and V - besides its dependence on the integration time of the 
system, Hence for given Tf , it can be assigned to each point 

of the R - V plane a value u (see 118) and consequently a corresponding 
degradation ratio ratio F(u) - F(R, V ) . Clearly then, the actual post- 
integration SNR will be a product of the expression for MSNR, given by 
(102) and the degradation ratio F(R, V t > , i.e., 


MSNR' = F (r, V t ) MSNR 


(130) 
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In the following we will examine the properties of F(R,V ). Comparing 

t 

(98) with (55), we find that the integration time of useful energy, ~X % 
which yields the desired MSNR as defined is given by 


/ 



v t * 


2Dr L 


R 







v t > 




2Dr n 


R 


(131) 


Figure 24 illustrates its distribution in the R-V fc plane. 



FIGURE 24. DISTRIBUTION OK THE T FOR 
FIXED INTEGRATION TIME. 
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The equation 



* constant (132) 


which assumes the form (127) , (131) 



Vt 

mk 


= c 


(for the shaded area in Figure 24) 


(133) 



(134) 


(for the region above the shaded area in Figure 24) 


determines in the R - V fc plane a set of curves with the property of 
connecting all points (R, V fc ) having the same degradation ratio* 

The curves for c ■ 0.5, c * 1.0 and c * 1.5 are shown in Figure 25. 
The value F(c) associated with each of them is the corresponding 
degradation ratio (compare with Figure 21) . Figure 26 gives the 
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FIGURE 25. CONTOURS OF THE DEGRADATION RATIO 


degradation ratio F(u) * F(R ^^R/2D T,) along and above the line 

v » A- R / 2D Tj. it can be shown that it oscilates around and for large 
R approaches the dashed curve in the Figure. 



FIGURE 26. RANGE DEPENDENCE OF THE DEGRADATION RATIO. 
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7*4*4 . Optimum integration titne of the system . As can be seen from 

Figure 25, the effect of the frequency deviation of the reflected signal 

2 

has no significance for ranges smaller than 8D l\ . Therefore, if the 

2 

required maximal range R « R is such that R < 8D l\ , then for 

& max max 7 

determining the optimal integration time and the corresponding 

transmitting power one can use the relations (109) and (ilC) . However 

2 

for maximal ranges greater than 8D the effect of the frequency 

deviation must be taken into account since part of the region of interest 
will enter in the area of high degradation if the integration time is 
chosen according to (109). Such a situation is depicted in Figure 27. 



FIGURE 27. AREAS OF HIGH DEGRADATION 
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If for instance, the maximal range is equal to 30D /X then because of 

the degradation ratio which for (R' V' ) assumes a value smaller than 
° max tmax 

0.1 (see Figure 26) it will be necessary to employ at least 10 times 
greater illuminating power than that calculated from (110) in order to 
achieve the indicated SNR. 

Before answering the question how we should determine the integration 

2 

time for ranges greater than 8D /X we will show that the MSNR' (106) 
below the dashed curve ABC in Figure 21 is not smaller than the MSNR' 
at the point C (or at the point A which has the same MSNR' as the point C) . 
Note that the point B is exactly at the crossing of the curves 
(see (133) f or u - C - 1 and R - R ' ) . 


v t = v tc = 


\JTxr; 


max 


t. 


v + - 


_ AR 




(135) 


*rotOT)UCIBlLlTY OF THE 
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The above will be proved if we show that (130) 


i\asnr’(r,vJi 


laloncj ABC 


r F/r,vJ MSNR^vJj 


|Qlonc| ABC 


^ v t mQ *) = o.8MSNR(R / mQX> V tc ) (136) 

For the part AB (.136) becomes 


F (R, Vj 


>. 0.8 


V* 




aloncf AB 


V, 


tmQx / 


oloruj AB 


(137) 


since the points between A and B lie above the line V * AR/2D 

u JL 

2 

where the MSNR is inversely proportional to and is independent of R. 
Introducing the expression for the curve AB, i.e., (Figure 27) 




max 


( 138 ) 
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in (113), we find that' we have to show 



(139) 


On the other hand (134) 



(140) 


But, (see Figure 21), o.8/u^ is equal to F(u) for u ■ 1, i.e. , at the 
point A and trailer than P(u) for any u > 1, i.e., for any (&, V) as we 

V 

are approaching th s point B. Hence (140) implies that (139) is 
satisfied and (136) is proved for the part of the curve between A and B. 
For the part BC (136) becomes 



( 141 ) 



i i-V- . , 

since the joists b et wee n B and C lie below die line V t " XR/2D ^ 

9 

/Where tie MB ie inversely proportional to £ and la independent of V . 
•/•'fib tie other land is (133) 



7 ? — 7 ===& * °* 6 

(h tfgsu ) 

U5r IT/ 


_R 

a* 


£ 


man 


(142) 


a 

Where again 0.8/a la equal to F(u) at the point C and aealler than 
F(u) for any other point between B and C. Thus (142) implies (141) 
and together with the previous result proves that the MSNR' on and 
below the curve ABC is equal to or greater than the MSNR* at the points 
A and C. 

flow, it is easy to see what will be the optimal integration time if 
we are to cover the area shaded w th vertical lines In Figure 28. 



FIGURE 28. R-V PLANE SITUATION FOR TH7. INTEGRA’ 1 T TIME WHICH 

MINIMIZES THE REQUIRED POWER TO COVER THE SHADED AREA. 
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As shown the integration tlae should be chosen is such s way that the 
contour u » 1 passes through the point ( > Mt i ▼ mnv ) > l.e., (see (133) 
for u " C • 1 and R » 1^) 

1 _ V tmo , <1M) 

This is because: (a) the Integration tine in 143 is the critical 

integration tine for the point (R » V ), that 1s, the one that 

maxi* ' * the USSR’ at that point (129) , and (b) the MSNR* associated 

with a^x toe points in the shaded area - and therefore with the 
points in the arer of interest - is greater than or equal to that at 
the point ( R saX » V taax ) 

To summarize, we give the general expressions for th *»ptinal 
integration tine and for the required illuminating power in terns of 
the maximal range, and the maximal tangential speed (at R * R ^ v ) , 

V 

tnax 


8D 2 - 


( xR m o, ? for R mw <^ 


tDV, 


t max 




( 144 ) 


]jz}^ R moy. 9 -f er R mQ y y 

v / A> 


tmw 



f 4(4») l (MswfU. MoRfDR«,\4 


wax 


0.8 Agjf 


* 4or 


^ ^t«»*n ” i 


^(WfMSNRl^R* V tmm f s6D » 

_____ ,ww x (1 , 5) 


lbs factor t.8 in the denominator of (145) is In cluded to account for the 

looses of the frequency deviation. It Is interesting to see that for ranges 

2 3/2 

greater than 80 /I the required power is proportional to & and is in- 

&8X 

dependent of the else of the array aperture D. 

Expression (145) can be written in the form 


*e({ P tmr> IT " ^ mox V tmox 9 Ur R„u* ^ 


(146) 


\ff v t m „ , for R„„ Z2g=ex 

*o A y 7 - 


(147) 


where 


is the resolution angle and 


_ 5(4Ttf(M5NR0 m ,„ N„R : 

*© — 

0, t (T 


(148) 


For a given angular resolution )( , wavelength A, gain of the 
transmitting antenna G^. and bistatlc cross-section of the targets, the 
equations (146) and (147) determine in the R-V t plane a set of curves wi'.h 
A ££p ^ as a parameter as shown in Figure 29. From these curves 




FIGURE 29. DEPENDENCE OF THE RANGE-TANGENTIAL SPEED COVERAGE 
ON THE INTEGRATION TIME AND THE TRANSMITTED POWER 
- EFFECTIVE RECEIVING AREA PRODUCT. 
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8. APPLICATION OF MICROWAVE IMAGING IN AN AIR-TRAFFIC SURVEILLANCE 
SYSTEM 

8.1 Sensitivity Characteristic of the Array Elements 

The cost of a Microwave imaging system with spaceborne illuminating 
source for the most part will be determined by the transmitted illuminat- 
ing power and the number of array elements. Therefore in order to save 
on these two components one should use array elements with gain as high 
as possible* or in other words, with a sensitivity characteristic shaped 
according to the solid angle of coverage. In the case of an air-traffic 
control and surveillance application, the solid angle of interest can be 
of the fora shown in Figure 30a. Consequently, dipoles with sensitivity 



FIGURE 30a. SOLID ANGLE OF FIGURE 30b. DIPOLE SENSITIVITY 

COVERAGE CHARACTERISTIC 

characteristic as that in Figure 30b can be considered suitable for build- 
ing the array. Since the angle a determines the gain of the elements 

m 

(fcr smaller o ffl the gain is bigger), its value will be a compromise be- 
tween the maximal gain and the maximal detectable elevation angle. In 
our discussions we will assume a ■ 20°. 


reproducibility of the 
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— = 2.^ smoC r 


(149) 


where: * gain of the element, Q «* solid angle of coverage, “ ef- 


fective area of the element, we find 


^ = -~i s A d = 


sw<rt 


2.TT 


(150) 


and for o ■ 20 


C» el =5.85 X 


Aj = -X 
* 2.15 


(151) 


8.2 Power Requirements 

From the discussions in the previous section, it can be concluded 
that the system can detect a target at any given range R < R^ a ^ provided 
that the transverse speed of the target is smaller than (146) (147) 


Arf(p ‘* ittH > 


‘f or ^max - ^ z 


*t max 


a o 2 .Ve \ I f O \ 8 X 

A e{(*t a — r — ~p/ z 5 f or > 2. 

*^o ^mQy 0 


(152) 


and the integration time is equal to (144) 


JL 5l!?j22L * /or R < - 8 ^ 

2 V ■» ‘'•non ^ - 




2=i ? I- M <153) 


't max 
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Our objective boro is to determine what the expressions (152) and (153) 
imply, in terms of transmitted power, number of array elements and duration 
of the integration time, under the circumstances of an air-traffic control 
ana surveillance application and for the following cases t 


Case 

I: 

Y " 

5 x 10 -3 

-3 

radians; 

X 

■ 0.2 meters 

Case 

XI* 

Y • 

2.5x10 J 

-3 

radians; 

X 

“0.2 meters 

Case 

III: 

Y " 

5 x 10 J 

-3 

radians; 

X 

“0.3 meters 

Case 

IV; 

Y * 

2.5x10 * 

radians; 

X 

“0.3 meters 


(154) 


2 __ —20 

Assuming: MSNR' ° 17 dB => 50, o * 10a , N = kTNF “ 10 Ws 

min o o 

R 53 4 x 10a, G ■ 250 Qq ( a 50-wavelength transmitting antenna dish) we 
find (148) 


k - 2.56 x 10 _6 Ws 
o 


(155) 


Expressing the effective area of the array as a sum of the effective 
areas of the individual elements, i.e. , (151) 

A eff - Si - "iTB* 2 < l5 « 


and introducing (155) and (156) in (152) one obtains 

► 

W < 


0.18169 -lO^N! 3 ££ -f » for < % 


0.18169 • 10 b NP t > f of <U > ~l 


t nte S 
K m<nt 


Letting 


(157) 


NP fc ■ q x 10 Watts x elements 


(158) 
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where q Is some dimensionless positive number, (157) becomes 


^tn>o< — - 
* 


0.18169 *»o" 

< 

0.18169 •10 M 

k 


Rmax ^ 

x^fax 


—■ 9 fcr ^ma, ^ x 2 
S 0 


(159) 


Combining (159) and (153) we rewrite the expression for the integration 
time in the fora 


( R <nox’ Himox) 


-10 p 2. 

0.e75l9*10 Jhnm. 

X x 


(160) 


Calculated values of V /q and qt_ for several maximal ranges, R , 

tmax jl ffiax 

and for each of the cases in (154) are given in Table 3. A graphical 



X»0.2; 

Case 1:^4 
Y* 

Y»5*10“ 3 
“ 64km ' 

X-C.i, y-2.5*10 3 
0% 

L?se II:^«256km 



X-0.3; y=5*10 3 
8X 

Case III;^2™96kni 

A«0.3; Y“2.5*10" 3 

03 j 

Case IV;^y3S4»fni j 

R 

max 

(km) 

V /q 

tmax 

(m/s) 

« T I 

(s) 

v. /q 

tmax 

(m/s) 

« T I 

(a) 

v„ /q 

tmax 

(m/s) 

H 

^tmax^ ** 
(m/s) 

m 

25 

145.360 

0.424 

72.680 

0.424 

327.060 

0.188 

163.530 

0.188 

50 


1.696 

36.340 

1.696 

163.530 

0.754 

81.765 

0.754 

75 

44.760 

3.816 

24.226 

3.816 

109.020 

1.696 

54.510 

1.696 


29.070 

6.784 

18.170 

6.784 

80.105 

3.015 

40.882 

3.015 


20.800 

10.630 

14.536 

10.600 

57.316 

4.711 

32.706 

4.711 


15.820 

15.264 

12.113 

15.264 

43.594 

6.784 

27.255 

6.784 


| 12.560 

20.776 

10.383 , 

j 

20.776 

i 

34.610 

9.234 

13.361 

9.234 


10.280 

27.136 

9.085 ' 

127.136 

28.328 

12.060 

j 20.441 

12.060 


TABLE 3. MAXIMAL DETECTABLE TRANSVERSE SPEED; REQUIRED INTEGRATION 
TIME DEPENDENCE ON TiiE DETECTION RANGE OF THE SYSTEM. 
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representation of (159) and (160), again for each of the cases in (154), 
la given in Figure 31. As can be understood, the descending curves cor- 
respond to V taax /<j, while the ascending curves correspond to qt^. The 

2 

dashed verticals denote the distance 8 X/y for the cases 1 and III. In 
order to illustrate how one can use the results in Table 3 or in Figure 31, 
we give an example. 

EXAMPLE: Suppose, an air-traffic control application requires detec- 

tion of aircraft according to the following specification: 



R 

V 

0 

- 25 km 

200 m/s 

25 

- 50 km 

150 m/s 

50 

- 100 km 

100 m/s 

100 

- 150 km 

50 m/s 


(Such a specification, where the required detectable transverse speed de- 
creases with the distance, is meaningful, since it is reasonable to expect 
that the angle between the aircraft path and the line of sight - for an 
aircraft approaching or leaving the airport - becomes smaller at greater 
distances.) From Table 3 we find for the first case (case I): 


Rv< 

25 km: 


145.36 ->q “ 

1.376 

R $ 

50 km: 

V cmax / ’ ' 150/q " 

72.68 “->q ” 

2.075 

P i 

100 km: 

V tm ax /<l ' 100/ ' 1 * 

29.07 ~>q - 

3.440 

R £ 

150 km: 

V t max 7 ’ ’ 50/<l ' 

15.82 =*>q » 

3.161 


(162) 


Consequently the required product of the number of elements and the trans- 
mitted power (134) will have to satisfy 

(NP ) , » q . x 10^ » 3.44 x 10^ Watts x elements (163) 

t min min 


REPRODUCIBILITY OF THE 
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(Por instance , If the transmitted power is 100 watts - the required num- 
ber of elements will be N ° 3440. In the *case of synchronous detection 


the required product NP t will be half as large; i.e. , ■ 1.72 x 10 J 

Watts x elements.) Going back to Table 3 with q ■ q ^ « 3.44*, we find 

that the maximal detectable speeds V _ , V , V , and V , _ ft 

tmax25 tmax50* tmaxlOO tmaxl50 

for the range domains (a) R4 25 km, (b) R 5 50 las, (c) R $ 100 km and 
(d) R < 150 km - with only one integration time in each case - and the 
corresponding integration times t i25 » t i5Q , t t100 , and t i150 are as fol- 
lows (see Figure 32): 


V t»ax25 " * 145 ‘ 36 * 500 " /s 
v tmax50 “ Vn * 72 -“ ' 250 «/• 
’'tmaxlOO ' U * 29 - 07 * 100 " /8 
’'taaxlSO ■U" 15 ' 82 * 54 - 4m/8 


for R $ 25 km 

for R ^ 50 km 

for R-5 100 km 

for R ^ 150 km (164) 


^25 “ 0.424/q <n ^ n * 0.123 seconds 
T I50 ° ^‘^^Snin " ^*493 seconds 
T I 100 “ 6,784/ Vn " 1,972 seconds 

Tjjjq ■ 15.264/q m ^ n • 4.437 second** (165) 


Comparing the diagram in Figure 32 (or the results (164)) with (161) we see 
that in order to cover the required region in the R - V plane we need 
only 3 different integration times, i.e., t hoO an< * T I150 ’ 

Repeating the same procedure for tue other cases we obtain the follow- 
ing results: 


The assumption is that *e tend to keep the produce Nl* t as small as possible; 
it is possible, of courru, to choose bigger q than q ^ . 
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Case II: (NP.) * * 5.504 x 10 5 watts x elements 

t min 

(2.752 x 10 5 - for synchronous detection) 

t X 50 “ °* 308 8 > V tmax50 “ 200 ^ 

T I100 " 1,233 8J V tmaxl00 “ 100 a/a 

t I 150 “ 2,733 85 V tmaxl50 " 66,7 m/s 

Case III: (NP^) ^ ■ 1.249 x 10 5 watts x elements 

(0.625 x 10 5 - for synchronous detection) 

t I 50 “ 0,604 8J V tmax50 “ 200 m/a 

T U00 " 2,414 85 V tmaxl00 " 100 m/a 

T I150 " 5,430 85 V tmaxl50 " 54,4 m/a 

5 

Case IV: (NP ) , « 2.446 x 10 watts x elements 

t'min 

C 

(1.223 x 10 - for synchronous detection) 

t I 50 ' 0,308 V tmax50 ' 200 Wa 
T H00 " 1,233 Bi V tmaxl00 “ 100 m/s 
t I 150 ' 2,773 si V t»axl50 ' 66,7 m/s 


8.3 Number and Position of the Beams 

The array is here assumed to be distributed over a horizontal 




(166) 


(167) 


( 16 b) 
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circular region. The radiation (or the sensitivity) pattern of a planar 
array focused at infinity is usually described in terms of the coordi- 
nates 8 and $ of a spherical-coordinate system, or in terms of the coordi- 
nates u and v of the so-called sinO -space. The relationship between the 
coordinates 6 and <J> of the spherical-coordinate system and the coordinates 
u and v of the sinG-space is given by 

u ■ rinG cos$ 

v ■ sind sin<j> (169) 

i.e., by the same relations which perform one to one imaging of the points 
on a unit hemisphere onto a plane (see Figure 33). 



FIGURE 33. PROJECTION OF POI. . ON A HEMISPHERE 
ONTO THE PLANE OF THE ARRAY 

A well known feature of the sinG -space is that the array pattern 
shape is invariant with respect to the direction of scan. This, and the 
reversibility of the transformation (169) which can be easily visuali7ed 
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by the representation in Figure 33 vill considerably simplify our task 
to determine the required number of beams and the angular position of the 
beams maxima; we can solve the whole problem by working only in the u - v 
plane where the closed curve connecting the 3 dB points of a given beam 
is invariant to the beam position. 

It can be easily shown that the 3 dB curve corresponding to a cir- 
cular array with equally distributed elements is a circle, of diameter 
Y “ X/D, in the u - v plane. Since we want the beams distributed in 
such a way that the distance between the neighbors is approximately the 
3 dB beamwldth, i.e., y; the problem of distribution of the beams reduces 
to problem of distribution of circles of diameter y. For convenience we 
will refer hereafter to the 3 dB circles as "beams". A simple distribu- 
tion of the beams - within the solid angle ft in Figure 30a - ie shown in 
Figure 34. A characteristic of this distribution is that while in the 
outside circle the neighboring beams just touch each other, the beams in 
the inside circle overlap. For o m ■ 20° (as was assumed befor ' the 
overlapping is nuwevei minimal since the inside circle is only 6% smaller 

-3 

(in diameter) than the outside. Dividing 2r by y ■ 5 x 10 and by 
-3 -3 

y ° 5 x 10 and by y * 2.5 x 10 wc find that one can associate 1257 
and 2514 azimuthal resolution cells with 200-wavelength and 400-wavelength 
arrays, respectively. Simila' Ly, dividing 1 - cosot^ " 0.06031 (for 

n » o _ o 

“ 20 ) by 5 x 10 and 2.5 x 10 we find that the region between the 
circles (1) and (2) in Figure 34 can be divided into 12 and 24 levation 
resolution cells for 200-wave] ength and 400-wavelength arrays, respectively. 
Thus, the total number of beams will be 
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200-wavelength arrays M * 1257 k 12 ® 15084 
400-wavelength array. M « 2514 x 24 ■ 60336 



FIGURE 34. DIS": cIBUTION OF THE BEAMS i'T THE SINO-SIACE 

Assigning to each beam a pair of indices (i,j) say in sue!* a way 
that the index "i" is associated with the elevation resolution cells and 
increases from the outside circle towards the inside circle, whereas the 
index ,, j H is associated with the azimuthal resolution cells and increases 
in the counter-clockwise direction starting from the upper side of the 
u-axis (for Instance the shao d bea*:, in the figure will be assigned 
indices (2,3)), one can express tne position of ihe maximum of any beam 
in terms of the coordinates w *» sinG and $ or in tct.iS of u and v. fas 
for a 200-wa r elcngth array one can write: 




( 171 ) 



— w- sin 





L * »... \TL 

j= 1... 1257 (172) 


TJve 8 as 2 expressions for a 400-wavelength array are 



i-0— k) 


a a.5* 10 


-5 


i = I-..Z4 



(173) 


Vi j=[l“( l 4) x2,<5 *l° 7 


i = l . . 4 2^ 

j = l.. . 2514 


(174) 


The elevation angles of the beams' maxima, for 2 00-wav el eng t 11 and 
400-wavelength arrays are given in Table 4. 

The 3 dB beamwidths - in elevation - at the smallest and at the lar- 
gest elevation angles are, respectively: 


200 wavelength array: 

5.732° (0.10004 radians) and 0.857° (0.01496 radians) 


400 wavelength array: 

4.052° (0.07072 radians) and 0.424° (0.00740 radians) 


( 175 ) 
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a (in degrees) 

Sj. (In degrees) 

4 

for a 200 wave- 

• for a 400 wave- 

X 

length array 

length array 

1 

4.052 

2.865 

2 

7.022 

4.964 

3 

9.069 

6.409 

4 

10.735 

7.585 

5 

12.177 

8.602 

6 

13.468 

9.512 

7 

14.647 

10.343 

8 

15.741 

11.113 

9 

16.764 

11.833 

10 

17.730 

12.512 

11 

18.649 

13.157 

12 

19.524 

13.772 

13 


14.362 

14 


14.923 

15 


15.474 

16 


16.002 

17 


16.514 

18 


17.011 

19 


17.494 

20 


17.964 

21 


18.423 

22 


18.871 

23 


19.309 

24 


19.737 


TABLE 4. ELEVATION BEAM CENTER LOCATIONS. 


8.4 Focusing 

From (87) and (94) above we conclude that the sensitivity of an end- 
fire focused circular array - along the direction of focusing - is given 
by 

5(R t » r f) - P sp R f ) 

Pr 



where; ■ array power outfjt, - focusing distance, * distance be- 
tween the radiating source and the array, 83 radiated power, ^ * constant 
of proportionality, \ ■ radiation wavelength and p Q ■ radius of the array 
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aperture. A graphical U lustration of the expression (176) is given in 
Figure 35. The curve (a) corresponds to Bj ■ while the curves (b) , 

(c), (d) and (e) correspond to Ip - Rg/2, 8^ - Rg/4, Rp * Rg/6 and Rp - Eg/ 8, 
where Eg la the byperfocal distance. As can be seen the focusing distances 
chosen in the Figure are the 3 dB distances. The sensitivity of the array 



FIGURE 35. ARRAY SENSITIVITY DEPENDENCE FOR ENDFIRE POINTING 
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to the hyper focal distance when focused at that distance Is taken as a 
normalizing factor (zero dB in the figure). 

The ordinates at R^/n, n ■ 1 . . .9, for the curves presented in 
figure 35, together with those corresponding to focusing at the hyper focal 
distance, are given in Table 5. 



V 9 

V 8 

V 7 

V 6 

V s 

V 4 

V 3 

V 2 

*H 

oo 

3.30 

3.05 

2.60 

2.24 

1.91 

1.44 

1.21 

0.659 

0.500 


3.86 

3.40 

3.06 

2.75 

2.25 

2.15 


B&a 

1.00 

V 2 

4.30 

4.0C 

3.82 

3.24 

3.36 


4.50 

4.00 

0.500 

V 4 

6.18 

5.77 

6.59 

5.93 

12.50 

16.00 

4.50 

0.659 

0.134 

V 6 

10.90 

10.55 

24.50 

36.00 

12.50 

2.64 

1.21 

0.360 

0.076 

V 8 

40.50 

64.00 

24.50 

5.94 

3.36 

1.44 

0.687 

0.250 

0.053 


TABLE 5. ARRAY SENSITIVITY VALDES FOR ENDFIRE POINTING. 

Since small elevation angles are of main interest for the air-traffic 
control application, i.e., almost endfire sensitivity is involved, we 
will examine the implications of the obtained results concerning such ap- 
plication. Going back to our examples we find (96) 


Case I: 
Case II: 
Case III: 
Case IV: 


Ry ■ 1865 meters 
Rjj ■ 7424 meters 
Rg ■ 2784 meters 
R„ * 11137 meters 


(177) 


We see that the distances within which the array has the ability of de- 
termining the range of the target by focusing are relatively small. If 
ranges smaller than those given in (177) are of little interest, then fo- 
cusing at distances other than infinity shouldn't be considered. If on 
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the other hand ranges saaller than those given in (177) are of interest 
and if a coarse range determination is of practical va 1 ue , then it will 
be useful to apply focusing. Thus if the minimal range of interest is 
1500 meters, then the useful focusing distances for the cases II and IV 
are: 


Case II: 

“fi - V 2 

S3 

3712 meters 


hi - V 4 

B 

1856 meters 

Case III: 

*n * V 2 

= 

5568 meters 


^2 " V 4 

= 

2788 meters 


“ra * V 6 

a 

1856 meters 


We repeat that the results above were developed on the assumption 
that the array detects from endfire direction - which of course is not 
completely true in detecting of aircraft. Larger disagreements are likely 
to appear only at smaller distances since only then do the elevation an- 
gles of interest become considerable. However substantial differences 
in the distribution of the 3 dB focusing distances cannot be expected 
since the hyperfocal distance by which this distribution is determined 
(97) is only slightly dependent on the elevation angle. 


9. SPATIAL CORRELATION EXPERIMENTAL EFFORT 

A modest experimental effort is proposed leading toward the verifi- 
cation of the microwave holographic array concept under development. 
Specifically, we propose to utilize the ATS-6 satellite after its return in 
Summer-Fall of 1976 tc the Western hemisphere together with suitable an- 
tenna-receiver modules to verify that an array of the type considered in 
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our Microwave Holographic Imaging System Design can be synchronized with 
the signal radiated from a synchronous-orbit satellite, permitting coher- 
ent combination of the signals received at the different elements of 
the array on the ground. 

The ground equipment will consist essentially of two receivers with 

* 

a conraon local oscillator, but separate 2m antennas. One receiving an- 
tenna will be fixed and one capable of being moved. The basic measure- 
ments taken will be with a steady sinusoidal L-band radiation from the 
satellite, measuring correlation of amplitude and phase shift variations 
of' the received signals. Measurements would be performed for a selection 
of receiver separations (vector position differences) . Once the signal 
strength, fading and phase coherence properties of the satellite-array 
geometry are established, the extent to which directional beams can be 
formed based on such synchronization can be assessed and limits on array 
size can be specified. 

The conduct of this experiment should fit in rather straightforwardly 
with the operational routine of the ATS-6 satellite. A convenient loca- 
tion of the receiving antennas will be at our Valley Forge Research Lab- 
oratory site in suburban Philadelphia (approximate location: 40°N. Lat., 

75.5° W. Long.). The laboratory is located at the top of a ridge (approx- 
imately 625 ft. elevation) with clear view in all directions. While the 
satellite was originally scheduled to be located at 95° W. Long, in the 
time period of interest (beginning Summer-Fall 1976), there is apparently 


A 2m antenna should provide about 28 dR SNR upon the satellite signal, 
assuming 4GW transmitter power and a 30 ft. antenna. 
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some doubt at present as to where it will be located. Line-cf-sight ele- 
vation angles for a variety of satellite longitudes have been calculated 
ae shown in Table 6 . For all these locations good visibility is achieved 
as shown. 


LONGITUDE 
(Degrees W. Long.) 

75.5 

90.5 

105.5 

120.5 
135.5. 


ELEVATION ANGLE 
(Degrees) 

43.75 

41.20 

34.40 

24.92 

14.10 


TABLE 6 . ELEVATION ANGLE FROM VFRC TO STATIONARY-ORBIT 
SATELLITE AS FUNCTION OF SATELLITE LONGITUDE. 


The experimental equipment to be assembled is illustrated by the 

block diagram of Figure 35 . The equipment is in general similar to the 

* 

radio camera array modules used for previous experiments at VFRC. The 
most important difference is the shift in operating frequency from 1068 MHz 
to 1550 MHz (the L-band down-link frequency of the ATS -6 communications 
subsystem). This accounts for the major component of the experimental 
equipment costs. 

As far as control of the satellite emission is concerned, all that 
ig required for an experimental run is to point the 30-foot dish at the 
Valley Forge Research Center site, then key on the L-band transmitter in 


* 


See VFRC Quarterly Progress Reports Nos. 1 to 15s also "Valley Forge 
Research Center Adaptive Array", B. D. Steinberg with E. N. Rowers, 
USNC/URSI Annual Meeting, Boulder, Colorado, October 1975, also "Design 
Approach for a High Resolution Microwave Imaging Radio Camera", B. D. 
Steinberg, Journal of the Franklin Institute, December 1973. 


REPRODUCIBILITY OF THE 
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FIGURE 36 . BLOCK DIAGRAM OF EXPERIMENTAL EQUIPMENT 

an unmodulated CW mode (either self-stabilized or locked to a stable 
transmission from one of the ground stations) for predetermined time in- 
tervals. It would appear sufficient to plan for total run durations of 
five or ten minutes during which the satellite would be keyed on and off 
at intervals of time ranging from one to ten seconds. Telephone contact 
should be sufficient between VFRC and the ATS Operations Control Center 
to alert operators as to readiness of equipment for new emissions. ("Rea- 
diness" consisting of moving one of the antennas, pointing it in the right 
direction, and checking st as of equipment and leads). 

Processing of data obtained should be very straightforward. Correla- 
tion characteristics of amplitude fluctuations with time and mean and var- 
iance of phase fluctuations are readily obtained. Careful measurement of 








- 99 - 


antenna spacing must also be carried out associated with the electrical 
amplitude and phase measurements made so as to permit assessment of ex- 
pected performance of a large array. Antenna spacings of the order of 
1000 feet in various directions can be provided conveniently at our Vi'RC 
site. 



